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MINIMUM WEIGHT OF FRAMES UNDER SHAKEDOWN LOADING 


Jacques Heyman! 
(Proc. Paper 1790) 


SUMMARY 


The first part of the paper discusses the minimum weight design of framed 
structures under both fixed and independently varying loads. A simple nu- 
merical example is given, for which the complete solution is obtained, and the 
results corresponding to the two types of loading are compared. 

An iterative method is then presented for the minimum weight design of a 


frame of any degree of complexity; this iterative method may lead to results 
which are slightly inexact. 


INTRODUCTION 


The type of frame to be considered in this paper is composed of prismatic 
members, joined together sufficiently rigidly so that the full plastic moment 
of any member may be developed at such joints. If distributed loads act on a 
frame, they will be replaced by a suitable combination of concentrated loads; 
under these conditions, plastic hinges can form only at joints and at loading 
points, owing to the piecewise linear distribution of bending moments between 
such critical sections. 

The normal assumptions of simple plastic theory will be made, of which 
the most important are that indefinite rotation can occur at a hinge position 
under constant moment (the full plastic moment), and that elastic and plastic 


deflections of the members are negligible compared with the overall di- 
mensions of the frame. 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1790 is 
part of the copyrighted Journal of the Engineering Mechanics Divizion, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 4, October, 1958. 

a. The results presented in this paper were obtained in the course of research 
sponsored by the Office of Naval Research under Contract Nonr-562(10) 
with Brown Univ. 

1. Fellow of Peterhouse, Univ. Lecturer in Eng., Univ. of Cambridge, Eng., 
Visiting Prof., Brown Univ., 1957-58, Providence, R. L. 
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The total weight of material used in a frame will be taken as proportional 
to the length and to the full plastic moment of each member; thus, if a typical 
member has full plastic moment My and length Mx, the total weight of materi- 
al in the frame will be given to some scale by 


ithe (1) 
where the summation extends over all the members. 

Minimum weight design under these conditions and with fixed loads acting 
on the frame was first discussed systematically by Foulkes;(1) the basic ideas 
can perhaps best be presented by means of a simple numerical example. 


Design Under Fixed Loads 


Consider the two span beam shown in Fig. 1, in which the point loads W, 
and Wp act at midspan. Each span has constant section, but the full plastic 
moments Mag and Mg may have different values. There are four possible 
modes of collapse of this structure, shown in Fig. 2. Modes (i) and (ii) occur 
for Ma < Mpg, and modes (ii*) and (iv) for Ma > Mg; the actual mode of col- 
lapse will depend on the relative magnitudes of W ae and Walp as well as 
upon the full plastic moments of the cross-sections. 

For example, if collapse occurs by mode (i), then, noting that the work done 
by the load Wa on a displacement of the collapse mechanisms is Wad Yay), 
while the work dissipated in the hinges is 3Ma 41, then 


(i) = 5 Wary (2) 


Similar expressions may be obtained from consideration of the other collapse 
mechanisms: 


(ii) My + 2My 

(iii) 2M, + My = (2 bis) 
(iv) 3M, = 


As a numerical example, suppose the two spans are equal, and that the 
loads have the magnitudes shown in Fig. 3, so that 


Wa%, = 1200 
(3) 
Wat, = 600 


(Note that it is not necessary to specify the units of length or of load in Eqs. 
(3)). Eqs. (2) then become 


ti) My = 200 
(ii) M, + 2M, = 300 
(iii) 2M, + M, = 600 (4) 
(iv) Mg = 100 
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Eqs. (4) are plotted in the design plane with axes Ma and Mp shown in Fig. 
4. According to the principles of simple plastic theory, a point in this plane 
which lies within the open convex polygon, defined by the mechanism lines, 
Eqs. (4), and shown shaded in Fig. 4, will represent a safe design of the 
structure, i.e. one capable of withstanding the given loads. A point on the 
boundary of the polygon represents a design of the structure which is just col- 
lapsing under the given loads; for this particular example, therefore, mode 
(ii) cannot occur for any values of Ma, and Mg. All points on the boundary of 
the polygon must now be examined to determine which represents the mini- 
mum weight design. 


From Eq. (1), the weight of the two-beam system is given by 


This weight line, when plotted in the design plane, is a line of 45° slope whose 
distance from the origin is proportional to the weight of the structure. Since 
the minimum weight is required, this will occur when the weight line is just 
“tangent” to the polygon. Thus the vertex M of the polygon in Fig. 4 repre- 
sents the minimum weight design for the two-beam system. It will be seen 
that two mechanisms ((iii) and (iv)) are operating simultaneously; from Eqs. 
(4), the equations of these two mechanisms lines can be written 


2aM, +a M, = 
(6) 
B Mp = 1008 


where multipliers a and 8 have been taken for the two equations. Adding Eqs. 
(6), 


2a M, + @+B)M, = 6000 + 100 B mn 


and comparing Eq. (7) with the weight Eq. (5), it will be seen that if a = 0.5, 
B =0.5, 


reo + M, = 350 (8) 

This superposition of alternative mechanisms is characteristic of minimum 
weight design; in the simple example above, positive coefficients a and 8 
could be found such that when these were multiplied into the mechanism 
equations, the equation resulting from their superposition gave the weight line. 
Interpreted geometrically, the condition that the coefficients a and 8 are posi- 
tive implies that the tangent weight line in Fig. 4 lies within the acute angle 
formed by the intersection of mechanism lines (iii) and (iv). This condition is 
necessary if the weight line is to be “tangent” to the convex polygon. 

These simple results were extended and proved rigorously by Foulkes for 
the case in which a frame is being designed to have N different full plastic 
moments. A geometrical representation can be made in an N-dimensional 
design space, the former mechanism lines being replaced by hyperplanes or 
“flats”. The mechanism flats for all conceivable modes of collapse define an 
open convex polyhedron (the permissible region). In general, the weight flat 
will be “tangent” at a vertex of the permissible region, each vertex being 
characterized by the intersection of, in general, N mechanism flats. Thus the 
minimum weight design will be characterized by N alternative mechanisms, 


. 
X = M, + 
‘A (5) 
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and these may be superposed with positive coefficients to give what may be 
called a Foulkes mechanism which coincides with the tangent weight flat. 

In exceptional cases, a vertex of the permissible region may give more 
than N alternative mechanisms; such behaviour does not interfere with the 
calculations. Alternatively, it may happen that one of the faces of the poly- 
hedron is itself parallel to the weight flat; for example, the weight line in Fig. 
4 might by accident have been parallel to mechanism line (iii), in which case 
a design represented by any point on the segment MN would have had the same 
minimum weight. (Mechanism (iii) in this case would itself have been a 
Foulkes mechanism; it will be seen readily that such behaviour will arise for 
the specific case of one span having twice the length of the other in Fig. 1). 

Returning to the more general example of Fig. 1, the four possible mecha- 
nisms of collapse shown in Fig. 2 may be superposed in pairs in an attempt to 
characterize all possible vertices of the permissible region. The significant 
possibilities are shown in Fig. 5, each of these mechanisms having two de- 
grees of freedom. Vertex M in Fig. 4, for example, is a particular case of 
mode D, corresponding to the superposition of mechanisms (iii) and (iv). For 
this mode D, from Eqs. (2), 


from which 


(10) 


" 

Ih 

= 
w 


Mp 
Now mode 4 can occur only for Ma > Mg; Eas. (10) give for this condition 
wt 
a+ >1 (11) 
Further, Eqs. (9) must be combinable with positive coefficients (a@ and § say) 
to give the Foulkes mechanism 


+4 gM, = (12) 


Thus 
2a 


A 
a + (13) 


from which 


a 


1 i 
B = zt + (14) 


The condition that 8 be positive gives 


2M, + = 

A Mp AA 

(9) 

— 

My 12 Wy B 

a 

4 

a= 

(15) 


A (i)+(ii) 


B (i) +(iii) 
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Inequalities (11) and (15) thus define the region D in Fig. 6 within which all 
minimum weight designs having Ma > Mg must lie. Similar considerations of 
the modes A, B and C in Fig. 5 lead to the correspondingly marked regions in 
Fig. 6. Fig. 6 is complete; knowing the ratios wal a/WB a and a/b, the 

minimum weight design can be determined immediately by using the figure. 


Shakedown Phenomena 


From Eqs. (6) and Fig. 4 for the simple beam example shown in Fig. 3, it 
is seen that the minimum weight design for this particular static loading con- 
dition is 
My = 250 


M, = 100 (16) 

Suppose now that the two loads Wa and Wg are not fixed, but may vary 
independently within the ranges 

O<¢ W, ¢ 1200 

(17) 

O<W, 600 
that is, the loads have the same respective maximum values as before, but 
can now take on any lesser positive values. The minimum weight design, Eqs. 
(16), is now inadequate to carry the system of varying loads, since increment- 
al collapse will for the corresponding load rapidly occur after a few cycles of 
loading. 

The phenomenon of incremental collapse may be explained briefly as 
follows: Suppose a given frame is acted upon by a set of loads varying be- 
tween prescribed limits, and that an elastic analysis is made for the effect of 
each load acting separately. For each critical section, the two combinations 
of loads can then be determined to give the greatest and least total elastic 
moment at the section considered. If the greatest and least elastic moments 
at all critical sections are less in absolute value than the corresponding full 
plastic moments, then the structure is safe. (“Full plastic moment” in the 
above sentence should really be replaced by “yield moment”, but this does not 
affect the argument). 

Suppose, however, that at one (or more) critical section the maximum 
elastic moment is greater than the full plastic moment. When the correspond- 
ing load combination is applied to the frame, yield will occur at that section, 
and a redistribution of bending moments will occur accompanied by the gener- 
ation of a system of residual moments (self-stressing moments) in the frame. 
Providing this residual system remains unaltered, the next time that the same 
loading combination is applied no further yield will occur. The response of 
the frame will be purely elastic, and the structure is said to have shaken down. 

However, plastic flow at some other critical section under some other load 
combination may have altered the residual moments; in this case, a second 
application of the original load combination may cause further plastic flow at 
the original section. Thus, as the loads vary independently, there is the possi- 
bility of plastic hinge rotations occurring at certain critical sections under 5 
one load combination, further hinge rotations occurring at other sections under 
another load combination, and so on. If this occurs, and if the sections which 
suffer plastic deformation are so located that they would form a plastic hinge 
system for a mechanism in the usual sense, then eventually the frame will 
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suffer severe deformations. Indeed, after only a few cycles of loading, the 
frame will look as if it had formed a mechanism under static loading. Thus 
failure by incremental collapse is in effect failure by the formation of a col- 
lapse mechanism; the hinges of this mechanism are not formed together, but 
some under one load combination, and some under another (or others). 


Minimum Weight Design Under Shakedown Loading 


Consider again the two span beam in Fig. 1, where now the loads Wa, and 
Wp can take any values in the ranges 


18) 
Wy < Pp (18) 


The moments of inertia of the two beams will be denoted I, and Ip; these 
quantities are obviously related to the full plastic moments Mag and Mg, and it 
will be assumed that 


I M 

A ‘An 

« 19 
M (19) 


A typical value for the exponent n might be 1.4; in most of the work which 
follows, n will be taken as unity for ease of computation, without seriously re- 
stricting the validity of the results. 

The following Table I gives the elastic bending moments at the four critical 
sections (1 to 4 in Fig. 1); the first two lines show the moments due to the 
loads acting separately, and these are then combined in the last two lines to 
give greatest and least values.Max and, Amin as the loads vary between the 
prescribed limits (18). Sagging bending moments are considered positive. 
The quantity k is the ratio of the stiffnesses of the two spans, i.e. 


4, 


(20) 


6 
32 | 


Table I 


Py ‘ J 
| 

“32 “Tne 32 “Tk 

3 Wptp 6 5+ 8k 

Me (8+ 5%) ° 

32 1+k 32 

32 32 32 1+k 
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Now, for any critical section i, the necessary and sufficient condition (2,3) 


for shakedown to occur is that a residual bending moment mj, can be found 
such that 


max 


min 21) 
My + ™ 2 


where My; is the full plastic moment of the member at section i. In the limit, 
the inequalities similar to (21) are replaced by equalities at a sufficient 
number of critical sections for an incremental collapse mechanism to be 
formed. 

Thus, suppose Fig. 2(i) gives the mode of incremental collapse of the two- 
span beam under the varying loads, hinges being formed at sections 1 and 2. 
Reference to Table I shows that 


(22) 


The residual moments m, and mg must be in equilibrium with zero external 
load; by virtual work, using the mechanism of Fig. 2(i) again, 


(23) 


or - 
Eqs. (22) and (23) solve to give 


(i) 


(24) 


and this equation must be satisfied if incremental collapse is just to occur by 
mode (i). Similar analyses may be made for the other three modes in Fig. 2. 
(It is not in fact necessary to introduce the residual moments since these are 


always eliminated. For example, Eq. (24) above could have been found by not- 
ing that since 


= 


max 


+m, = -M, 
then the sum obtained by adding the first of these equations multiplied by 
(264) to the second multiplied by (- 6) will give precisely Eq. (24), m, and 
mg disappearing by virtue of Eq. (23). The quantities (261) and (- 6;) are the 
hinge rotations of the mechanism considered, Fig. 2(i)). The other three 
modes give 


(iii) 


(24 bis) 


Me + ™ 
+m = M 
32 A 
32 32 1l+k 2 A 
(20, )m, + (-0,)m, = 0 
= A + q 
| 
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For the simple numerical example, Ya = Ip = 1, suppose again that 


October, 1958 


P = 1200 


Then Eqs. (24) give 
(i) 


(ii) 


(iii) 


= + 100 
(iv) Mp + 


where k, for the particular case 4, = ip, is given from Eqs. (20) and (19): 


(27) 


For n = 1, Eqs. (26) are plotted in Fig. 7, in a way similar to that for the 
static loading, Fig. 4. Comparison of the two figures shows the general re- 
semblance of the cases of fixed and of varible loads; minimum weight in Fig. 


7 is again represented by the point M at which two alternative mechanisms 
form, and for which 
My = 277 


= 123 
(28) 


Note the marked increase in weight X, from 350 in Eq. (8) to 400 for the varia- 
ble loads having the same maxima as the fixed loads. 

The lines giving the boundary of the permissible region in Fig. 7 are slight- 
ly curved. There is thus the possibility of the weight line being tangent to the 
permissible region at a point which is not an apex. Fig. 8 illustrates this. 

For the particular ratios aN = 1.621 ip, and Pala = 2Pp B, the straight line 
joining the apices MN of the permissible region is parallel to the weight line. 
Since the mechanism line (iii) is slightly curved, the true minimum weight 
will occur at a point on the boundary of the permissible region between M and 
N. The weight of the designs corresponding to points M and N is 561.7 (for 
Pala = 1200, as before); the true minimum weight is given by calculation as 
561.1. The curvature of the portion MN of mechanism line (iii) in Fig. 8 is so 
slight that the point of tangency cannot be distinguished in this figure. Fig. 9 
shows the calculated weight X for designs lying on the line MN, and gives the 
minimum weight for the ratio Mp/Ma equal to about 0.75. Thus confining at- 
tention to the apices of the permissible region leads to an error of about 0.1%. 
In the following, no attempt is made to allow for this slight error, and it is 
assumed that minimum weight always occurs at an apex of the permissible 


region. 


an 
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A brief discussion of the full solution for the two-span beam problem will 
now be given. As before, Fig. 5 gives the significant possibilities of alterna- 
tive collapse mechanisms corresponding to apices of the permissible region. 
For mode D of incremental collapse, for example, (combination of mechanisms 
(iii) and (iv) as in Fig. 7), Eqs. (24) give 


be 96 1+k ‘ 


where k is given by Eqs. (20) and (19), 


(30) 


As before, mode D can occur only for Ma 2 Mg; from Eqs. (29), this condition 
gives 

> 5+8k 
85k (31) 


In the limit when Ma = Mpg, then k = Va/ip from Eq. (30), and the equality 
(31) has been plotted in Fig. 10 dividing the regions B and D from A and C. 
The “vertical” dividing lines in this figure are more difficult to determine, but 
implicit analytical expressions may be obtained; the lines shown in Fig. 10 
are drawn for the exponent n = 1, as usual. Figs. 6 and 10 again bear marked 
resemblances, the straight lines being replaced by slightly curved lines. 


General Iterative Method for Shakedown Design 


It will have been noted that shake-down design is essentially an elastic 
problem, since the elastic response of the structure must be used. Now the 
elastic bending moments in the structure depend on the stiffnesses of the vari- 
ous members; unfortunately, for the design problem, these stiffnesses are not 
known at the start of the calculations. Thus the vicious circle obtains in which 
it is necessary to know the sizes of the members in order to obtain the elastic 
response, and to know the elastic response in order to determine the sizes of 
the members for minimum weight. It is true that by employing some device 
such as a stiffness matrix (written in symbols), a complete formal solution 
can be obtained; such a solution is unwieldy, however, for all but the very 
simplest frames. 

The iterative method proposed assumes a certain design of the frame. The 
method of moment distribution, or some other convenient means of elastic 
analysis, will then enable the elastic response under the given loads to be de- 
termined. In particular,.4 ™€X and.Mmin can be calculated for each section. 
These values ™ax and.4{min are then taken to be unvarying, and the 
minimum weight design is derived under these conditions. This minimum 
weight design will differ, in general, from that originally assumed. New 
values of:M max and. Mmin are then calculated, and the process is repeated, 
convergence usually being reasonably fast. 
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The computations will be made in the tabular form proposed recently. (4,5) 

As an example, the same two-span beam problem will be considered, with 

A= le = 1, Paha = 2Pp kp = 1200, k = (Mp/Ma)®" (see Fig. 7). The compu- 
tation will be started by assuming that one beam is of negligible strength com- 
pared with the other, say Mp << M4, i.e. k = 0. The elastic bending moments 
may then be computed directly, or, for this simple example, k = 0 can be 
substituted in the expressions for 44MaX andM Min in Table I. The resulting 
values of & max and M min have been entered as Step O in line 4 of Table II, 
the first four values corresponding to./A-M€X at the sections 1 to 4. 

The three preceding lines in Table II, labelled “data”, represent the hinge 
rotations of independent elementary mechanisms from which all conceivable 
mechanisms can be obtained by linear combination. For example, line 1 in 
Table II represents the mechanism in Fig. 2(i); line 10 in the same Table, 
which is the sum of lines 1 and 3, represents the mechanism in Fig. 2(iii). 
These data equations also give the coefficients of the independent homogeneous 
equilibrium equations which must be satisfied by any system of residual 
moments; for example, if residual moments m, and mg act at sections 1 and 
2, then line 1 of the Table states that 


(32) 


which is identical with Eq. (23) above. 

Examining line 4 of Table II, and remembering that it is assumed for the 
moment that this line represents a fixed elastic response, it will be seen that 
a beam system having 


My = 300, M,= 112.5, X = 412.5 (33) 
is safe. (The largest absolute values of the elastic moments in spans A and 
B are asterisked in line 4). The most general system of residual moments 
will now be constructed in an attempt to modify line 4 in such a way that the 
total weight X is reduced, and subject to the following condition: If the maxi- 
mum moment in span K is Mx, then the sections at which Mx occurs will be 
“attacked” by residual moments -mx. Thus all “hinges” (sections at which 
the moment My, acts) will be maintained while the value of Mx is reduced in 
an attempt to reduce the total weight of the structure. 

For the numerical example, it will be seen that the largest moment in span 
A occurs at section 1. Thus the residual moment m, will be set equal to -m,; 
similarly, the largest moment in span B occurs at section 3, and mg will be 
set equal to +mp. Data Eq. (3) requires that Mg = Mg = Mp, while data Eq. 


(2) gives mg = ; m3 = 3 my. These values are displayed in line 5 of Table I. 


Further, data Eq. (1) requires 2m, + Mp = 0, so that line 6 of Table II gives 
the most general distribution of residual bending moments that can be used to 
modify line 4. For example, the moment at section 1 can be reduced to 

300 - ma, representing a saving of material of 1 x ma, while the moment at 
section 3 is increased to |112.5 + 2ma]|, representing an increase in material 
consumption of 1 x 2ma. Such a change implies a net increase in material 
consumption; evidently if ma, is taken negative, a saving in material will be 
made. The maximum saving occurs for ma = -6.2, noted in line 6 of Table II; 
line 7 gives the residual moment distribution to be added to line 4, and line 8 
gives the result of this operation. m, was limited to the value -6.2 since a 
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new hinge is formed at section 4 for this value. No further change is possible 
in an attempt to save material, and Step I, occupying lines 5 to 8 in Table II, 
is terminated. (In general, Step I may be repeated by making suitable modifi- 
cations; an example is given below). 

Step II, lines 9 to 13, checks that the Foulkes mechanism which always ap- 
pears at the end of Step 1 is a true Foulkes mechanism. Line 9 gives the signs 
of the hinge rotations corresponding to line 8; note that the negative hinge has 
been transferred from the second half of the Table. The original data 
equations are now combined in such a way that zero coefficients occur at criti- 
cal sections (i.e. section 2 in this example) where no hinge has been formed. 
This has been done in lines 10 and 11; line 12 gives the most general combi- 
nation of the two mechanisms under coefficients a and 8. For a true, as op- 
posed to a false, Foulkes mechanism, the quantities in brackets in line 12 
must be positive, and, in each span, must sum to the length of the span. Line 
13 displays the hinge rotations under these conditions, and it will be seen that 
the solution is satisfactory. (Had one or more of the quantities in brackets 
been negative, corresponding to a false Foulkes mechanism, then Step I could 
be re-initiated by dropping the requirement of maintaining the hinge at one of 
these false hinge positions.) 

From line 8 of Table II it will be seen that subject still to the (false) as- 
sumption that the elastic response is invariant and given by the bending 
moments in line 4, a safe design of the two beams is 


My 306.2, Ma 100.0, X = 406.2 (34) 
These values give k = 100.0/306.2 = 0.326, and this revised value of k will give 
rise to a modified set of elastic moments in the frame. This new elastic re- 
sponse is entered as a new Step O in line 14 of Table IJ; as before, the total 
weight of the beam system is reduced, and the design corresponding to line 16 
is 


M, = 283.4, = 118.6, X = 402.0 


The new value of k is 0.418, and the process is repeated; after four iterations, 
the values 


My = 277 4, Mp = 122.8, xX = 400,2 


(36) 


are derived, compared with the correct values given in Eqs. (28). 

If k had been taken as (Mp/Ma)®", where n is not equal to unity, it is evi- 
dent that no complications would have arisen. The iteration may have been 
shorter or longer, depending on the starting value of k. 

A second example will show that errors can be introduced into the analysis 
by using the concept that the elastic response is fixed for any one iteration. 
Consider the same beam system, Fig. 1, but with 


P,t, = 1200 
4, = 1.8 
4. = 1 


Reference to Fig. 10 shows that the design point lies within the area B, for 
which the two beams have equal sections, and exact calculation gives 


(37) 
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M = M = 21302, » 4 = 597.0 


A B (38) 


Now k = 1.8 Mp/Ma: the iterative process will be started assuming a uniform 
beam, i.e. k = 1.8. The working is given in Table IlI(a), and after three iter- 
ations line 10 gives 


M, = 260, MB = 137, X= 605 (39) 
This design is satisfactory in that the corresponding beam system will just 
fail by incremental collapse under the given loads; however, comparison with ' 
Eqs. (38) shows that it is not the minimum weight design, the weight being just ‘a 
over 1% too high. q 

The iterative method has converged on the wrong vertex of the permissible 
region; using Fig. 7 for illustrative purposes (this figure is drawn for 

A= Ye and not Qa = 1.8 Lp), the iterative method has given vertex M as the 
minimum weight design whereas the true minimum occurs at vertex N. The 
error (which, as was seen, was small) has been introduced because of the 
false assumption that the elastic response may be regarded as fixed while the 
ratio Ma/Mg is changed. Line 2 of Table III shows the saving consequent on 
changing ma by a certain amount; if mg is given the value unity, for example, 
a saving of 1.8 units of material will be effected in span A, but an increase of 
2 units of material will occur in span B. Thus Ma was chosen negative to 
give a net saving of 0.2 times the value of mag. In effect, the iterative method 
assumes that the slope of mechanism line (iii) in Fig. 7 is constant and equal 
to 2, that is, greater than 1.8; in fact, the mechanism line is curved, and in 
the region MN has a slope less than 1.8. 

Table II(b) shows the result of choosing ma positive; line 1 has been 
transferred to line 11, and line 13 shows Ma = Mg, thus concluding the iter- 
ation. Between lines 11 and 13 the weight has increased from 594 to 597, but 
the latter value is less than that given in line 10. The iterative method will 
tend to break down in this way when the slope of the Foulkes mechanism line 
(i.e. of the weight line) is close to the slope of one of the lines representing 
an actual mechanism of failure. Thus a design point in Fig. 10 which lies 
close to the lines separating region B from D, or A from C, will represent a 
case for which the iterative method may lead to errors. (The iterative method 
will replace these curved dividing lines by straight lines at Ma/dp = 2 and 
1/2 respectively). The danger signal is given, as in Table III, by a large 
change in one of the attacking moments (ma) producing only a small change in 
total weight of the structure. If this happens, then the designer should check, 
as in Table III(b), whether an apparent increase in weight in fact leads to the 
correct solution. 


Final Example 


As a final example, consider the simple portal frame shown in Fig. 11 in 
which the two columns have equal sections, possibly different from that of the 


beam. If the loads shown are static, then the minimum weight design occurs 
for 


X = 30 My + 20 M, = 3000 
M, = Mp = 60 (40) 
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For the shakedown design, the two loads will be assumed to vary inde- 


pendently between zero and the maximum values shown. The stiffness ratio 
constant will be taken as 


I Mp 1,4 
k= = (41) 
A A 
The iteration will be started with a uniform frame, i.e. k = 3/4. Table IV 
gives the working. This table is split into two haives, as usual, the left hand 
half dealing with the maximum and the right hand half with the minimum 
elastic moments, 

The four independent mechanisms of collapse are displayed as line 1 to 4 
of the Table; line 1 gives the sidesway and line 2 the beam mechanism, while 
lines 3 and 4 correspond to joint rotations at the two corners of the frame. 
Line 5, Step O, gives the maximum and minimum elastic moments for a frame 
with equal beam and columns (i.e. k = 3/4). The two asterisked moments are 
the largest which occur in the respective spans for either half of the Table. 
ot In line 6 is constructed the most general system of residual moments which 
ag . can be used to modify the moments in line 5. A moment +m, “attacks” 

section 3, where the largest negative moment in span A occurs, and -mp at- 
tacks section 6 in span B. Since the frame has three redundancies, a supple- 
mentary moment x is required to complete the residual distribution. Lines 
7, 8 and 9 form an operations table; each line may be used independently of 
the others. Thus a unit change in my will save 30 units of material, while a 
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unit change in mp will save 20 units. If mp is taken as 5.45, line 10, and this 
line is superposed on line 5, line 11 results, which shows that a new hinge is 


WEIGHT OF FRAMES 1790-23 


formed at section 4. The total material saving is 109 units, which has been 
noted in line 8. Had line 7 been used, this also would have resulted in the 
same saving of 109 units. 

Step I requires that all hinges formed during the operation should be main- 


tained. Thus the newly formed hinge at section 4 must be attacked by -mag. 


Now the total residual moment at section 4 is (2m, + 2Mp + x); thus 


em, + (42) 
Eq. (42) enables one of the residual moments to be eliminated from the oper- 
ations table; lines 12 and 13 show the result of eliminating the supplementary 
moment x. Step I can now proceed; use of line 12 saves 206 units of material, 
and hence is more efficient than use of line 13, which saves only 56 units. A 
new hinge is formed in line 15, mg is eliminated in line 16, and finally no 
further move can be made by Step I at line 18. 

Lines 19 to 23 check that a true Foulkes mechanism has been formed, and 
it will be seen that all the hinge rotations in line 23 are of the right sign. As 
mentioned earlier, had one of the numbers in brackets been negative, Step I 
could have been reinitiated by dropping the requirement that this false hinge 
be maintained. 

Line 18 gives the following minimum weight design, on the assumption that 
k = 3/4: 


M, = 58.7 = 73.6 X = 3233 
From values (43), however, k from Eq. (41) is found to have the value 1.03. 
Remembering the general behaviour exhibited in the example worked in Table 
II, the next iteration is made for k = 1.1. Lines 24 to 36 of Table IV summa- 
rize the calculations, and lead to 


My = 5607, Mp = 7603; } 4 = 3227, k = 1.13 


(44) 
Evidently the iterations are converging to the final design: 


My, = 56.6, Mp = 76.4%, X= 3226, k = 1.1% (45) 
No errors have been introduced in this particular solution by the non-linearity 
of the basic equations. The minimum weight design under shakedown con- 


ditions uses 7-1/2% more material than the corresponding design under static 
loads. 
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Proceedings of the American Society of Civil Engineers i 
EFFECT OF END-FIXITY ON THE VIBRATION OF RODS q 
D. Burgreen,! 
(Proc. Paper 1791) 
A study is made of the free vibration of rods with end supports of equal i 
elasticity. Expressions are derived relating the end-fixity to frequency and _ 
mode of vibration. With the proper type of negative springs attached to the es 
ends of the rod, it is shown that vibrations can be completely surpressed. , 
The relationship of this phenomenon to the instability of beams with nega- TR 
tive end-fixity is discussed. a 
NOMENCLATURE 
Young’s modulus of elasticity 3s 
moment of enertia of bar 
p mass density of bar Pe 
A = cross-sectional area of bar = 
X, = f(x), = equation of mode 4 
W,, = natural frequency of vibration a 
Bn = (pAw,?/ED! 
L = length of bar a 
a@ = proportionality constant defining end fixity used in the expression 


depth or thickness of bar 


maximum amplitude of vibration, or maximum deflection 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1791 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 4, October, 1958. 
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€, = maximum extreme fibre strain 
k. = amplitude-strain factor defined by : =k, w/t 


w  =distributed beam load 


In a vibration program dealing with the effects of high velocity water 
flowing parallel to long thin rods, one of the problems that arose and re- 
quired investigation was the effect of the end connections of a prismatic test 
rod on the natural frequencies and modes of vibration. Specifically, the 
questions that required an answer were: 

1. If the fundamental frequency of vibration of the rod with an equal 
amount of end-fixity at each end, is known, what are the higher natur- 
al frequencies of vibration? 

2. If the extreme fibre strain in the rod is known, what is the amplitude 
of vibration? 

Texts on vibration discuss the problem of vibration of prismatic bars 
with specific end conditions such as pin ends, built in ends, or free ends, 
and generally list some of the lower natural frequencies of vibration of these 
bars. In practice, however, most bars or beams are neither pin-ended nor 
built in, but have end connections with some amount of elasticity in them. 

If the amount of elasticity in the end connections is known, then the natural 
frequencies and modes of vibration can be determined. 

It will be shown that a knowledge of the amount of end-fixity in the rod 
supports is not necessary in order to answer the above questions. The end- 
fixity is simply the parameter which relates the desired quantities. There 
are, however, occasions when the end-fixity as determined by a measure- 
ment of the fundamental frequency of vibration or from strain deflection 
data may be of direct use. For example, it may be desired to determine 
accurately the load carrying capacity of a beam in a structure—which de- 
pends upon the amount of end-fixity of the member. This can be done by the 
use of the curves in Fig. 2 wherein the end-fixity is plotted against the 
natural frequency of vibration. Some characteristics of elastically restrain- 
ed beams that are of practical interest are discussed in the text that follows. 

The equation for the vibration of the prismatic bar with elastically built 
in ends, shown in Fig. 1 is 


+4 pA O%y = 0. 


(1) 
at? 
The s»lution to this equation is taken as 
y= > Xx, sin wt, (2) 
n=1 
which gives the equation for the mode of vibration as 
2 
4009 
x, =0. (3) 
With the notation EI 
4 2 
= PAwn , 


EI 
Equation (3) becomes 


eek 


Xn - pa xX, =0. (4) 
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The solution to Equation 4 is 
X=Asin Bx+Bcos Bx+C sinh Bx + D coshBx (5) 


The ends of the bar are assumed to be elastically built in with the angle 
of rotation proportional to the applied moment. This is the common type of 
linear torsional elasticity expressed as 


M = constant 6. 
The boundary conditions representing this are expressed as 
Big yg (6) 
and 
ad 
The other boundary conditions are 


a 7 
Xi (7) 


Setting Boundary Conditions (6) and (7) into Equation (5), the following 
frequency equations are obtained: 


tan BL/2 + tanh BL/2 (8) 
aL = 2pL 


cot BL/2 - coth BL/2 (9) 


The above expressions refer respectively to the symmetrical and to the 
anti-symmetrical modes of vibration and are plotted in Fig. 2. These curves 
give the relationship between the end-fixity and the natural frequencies of 


‘ 
vibration. The curve designated as n=1 refers to the first symmetrical mode i 
of vibration of a prismatic bar whose end-fixity is a@. The special case of a 
pin ended beam with zero end-fixity is included here as the point of inter- ae 
section of the curve with the BL axis. At this point aL is zero and BL=7, a 


which corresponds to the known fundamental frequency of vibration of a pin 
ended beam. The higher natural frequencies of vibration of a pin ended beam 
are given by the intersection of the § L axis with the successive curves. 

The upper part of curve n=1 approaches an asymptote at 8 L=4.73; that 
is, as aL approaches infinity, signifying a bar with infinite end restraint or 
built in ends, BL, the root of the Frequency Equation (8), gives the fre- 
quency of vibration of a built in beam. The asymptotes, as aL approaches 
infinity, on the higher numbered curves, similarly give the higher frequen- 
cies of vibration of prismatic bars with built in ends. 

On curve n=1, which refers to the first symmetrical mode of vibration, 
it may be noted that low order wave shapes are included up to that of a beam 
with built in ends (one full wave). Curve n=3 which represents the second 
symmetrical mode of vibration starts at al7*- © with the mode shape of a 
a - beam with built in ends, or one full wave, gradually changes its shape with 


increasing BL to become a pin ended beam in its second symmetrical mode : se 
with one and a half waves, at a L=O, and finally at @L~ oo assumes the shape Li 
of the second symmetrical mode of a beam with built in ends, or a beam ian 
having two full waves. The progress of change in shape with end-fixity and a” 
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frequency may be traced in the same manner for the even numbered anti- 
symmetrical modes. 

Of particular interest is the significance of the negative values of the end- 
fixity. This is the case in which the springs at the extremeties of the bar, 
instead of offering restraint and acting as a restoring force against the 
inertia forces of the vibrating bar, actually assists the inertia forces and 
act in opposition to the restoring forces due to the bending rigidity of the 
bar. 

A bar with built in extremities offers the greatest end restraint. As the 
end-fixity is gradually relaxed by attaching the ends to successively softer 
springs, the restoring forces at the bar ends is diminished until the limit- 
ing case of a bar with zero end restraint, or one with pin ends, is attained. - 
Any lesser end restraint must now be a negative end restraint. It may be 
possible to attach devices to the extremities of the bar which can be called 
negative springs, which cause the bar to have greater end slopes, in the 
absolute sense, than it would have if it were pin ended. Such devices assist 
the motion when the bar is moving from its straight line position and hold 
it back when it is moving toward its straight line position, and therefore 
have the effect of reducing the frequency of vibration. 

A method of attaining a negative end-fixity is shown in Fig. 3. Steel 
pieces are attached to the ends of the rod, supported by pins, as shown. For 
the purposes of the illustration, the mass and moment of inertia of the steel 
attachments must be disregarded. Two stationary magnets are places in 
close proximity to the end fittings, so that the bar in its straight line po- 
sition, as shown in the figure, will not have any tendency to deflect either 
upward or downward. It can be seen in this arrangement that when the bar 
deflects downward, the effect of the magnets is to increase this deflection, 
and when the bar deflects upward, the effect of the nagnets is to increase pis 
the upward deflection still further. The magnets also restrain the movement 
toward the straight line position. It is seen then that this behavior is exact- 
ly in opposition to the behavior of normal springs. It is apparent that when 
strong magnets are used a small deflection from the straight position may 
cause the bar to buckle. That is, there may be a critical magnet strength 
associated with instability of the bar in Fig. 3. 

When BL is set equal to zero in Equations (8) and (9), signifying zero 
frequency of vibration, @L=-2 and a@L=~-6 are found to be the end-fixities 
associated with zero frequency of vibration in the symmetrical and anti- 
symmetrical modes of vibration. It is well known that the conditions for 
zero frequency of vibration are often the same as the conditions required 
for instability. This phenomenon is usually associated with axial loads 
rather than negative springs. It is of interest to determine whether this 
correspondence exists also for negative springs. 

Consider the bar in Fig. 4, which is the same as the bar under discussion, 


except that now it is uniformly loaded. Using the same boundary conditions 
as before, the equation for the deflected bar is 


wxt = wL2x2 ( 6 | +aL ) (10) 
24 48 2+ aL 384 2 +aL 

Note in the above equation that when a =O, the deflection equation for a pin 
ended bar is obtained, and when @~+.o, the deflection equation for a built in 
bar is obtained. The particular feature in the above equation that is of in- 
terest is the indication of instability when aL approaches ~2. This is the 
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critical value of the magnet strength that would produce an instability con- 
dition. In Fig. 5, the same bar is shown with an antisymmetrical loading 
expressed as 


W-=wsin 27x 


L 
: and the resulting deflection equation is a 
Ely - WL gin 2H cw ) ai) off 
6+aL 16m \ 64 aL 
The critical vaule of the end-fixity in this equation is aL=~6, indicating " 
i that this value of end-fixity refers to an anti-symmetrical type of instability. * 


The correspondence of instability due to negative spring supports, of end- 
fixity a@L=-6 when a bar is antisymmetrically loaded, and of end-fixity of 
otL="2 when the bar is symmetrically loaded, to zero frequency is thus 
shown. 


w sin 


L 


This may also be seen from the n=1 curve in Fig. 2. When the end-fixity 
of a bar vibrating in its first symmetrical mode reaches the limit a@L=-2, 
then no vibration is possible. The second curve, n=2, in Fig. 2, shows that 
the natural frequency of vibration in the first antisymmetrical mode ap- 
proaches zero as the end-fixity aL approaches the instability value ~6. The 
higher numbered curves corresponding to higher mode shapes do not have 
any limiting values of negative end-fixity which exclude the possibility of 
vibration. The negative values of end-fixity of the higher modes have, how- 
ever, limited physical significance, since it is not generally possible to get 
past the lowest end-fixity that causes instability, that is, aL=~2. 

Of practical interest in the study of bars with elastically supported ends 
is the relationship between amplitude of vibration and strain. The constant 


k, which relates amplitude and extreme fibre strain is given as . 
k 6 t 


k may also be expressed as, k=2K/L2x” . In experimental tests, data is 
sometimes recorded by means of strain gages in terms of extreme fibre 
i ‘ strain, which must be converted to amplitude. The above relationship which 

is valid for loaded bars as well as for vibrating bars is useful in experimen- 
tal work. If the value of k at the center of the bar is defined as kK, , then for 
a uniformly loaded beam with elastic end supports it is given as° 

_ aL+ 10 (12) 
ou =68aL + 48 
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and for a beam with elastic end supports having a concentrated load at the 
center of the span it becomes 


k = aL + 8 (13) 
120L + 48 


the above expression may be used, for example, in determing the amplitude- 
strain relationships of a uniformly loaded beam with undetermined end- oF 
fixity when calibrations are made with a concentrated load. For a bar which Se 
is vibrating freely, a similar expression may be derived. The equation for 
the symmetrical mode has the form 


X =D (cosh Bx - cos Bx) (14) 


which satisfies the condition of zero displacement at the ends. For the 
center of the beam this gives 


cosh BL/2 
| k | Xy-9 = COS gL/2 (15) 
cos BL/2 


and relates the amplitude of vibration with the extreme fibre strain at the 
center of the beam as given below 


cosh BL/2_ 
6 cosh BL/2 L 


€ 2(BL/2)" 14 coshpL/2 | (16) 
cos B 2 


The above expression for k. does not include the end-fixity @ directly 
as in the case of beams. The deflection-strain relationship obviously de- 
pends upon the rigidity of the end supports and its dependence upon @ is 
seen in Equation (8) wherein BL is given as an implicit function of a L. 
Equation (15) is plotted in Fig. 6. It may be observed that the ratio of de- 
flection to extreme fibre strain becomes smaller as the frequency of vi- 
bration (or B L) increases. It signifies the crowding of a large number of 
waves into the rod length, which must necessarily produce high curvatures 
and small deflections. A limiting case is represented by a pin ended beam 
which vibrates in a half sine wave. The amplitude-strain factor for a sine 
wave is 2/12. Fig. 6 gives BL=7 for this value of k,, and it corresponds 
to the known fundamental frequency of vibration of a pin ended beam. The 
curve also shows that k, approaches a limit of 1/4 as the frequency ap- 
proaches zero. It was shown that the natural frequency of vibration of the 
first symmetrical mode of a bar approaches zero when the end-fixity 
approaches its instability value of @L="2. Since no vibration occurs when 
aL=2, the bar may simply be considered uniformly loaded by its own weight. 
Equation (12) which shows the variation of kK, with the end-fixity a, for a fs 
uniformly loaded beam verifies that the value of the amplitude-strain factor ‘@ 
becomes 1/4, as aL approaches ~2. In Fig. 6, it is seen that 1/4 is the 
largest value that ky may attain, indicating that the shape in which the bar 
becomes unstable is that of least curvature at its center. 
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In the practical application of strain to amplitude conversion it will be 
found that the strain-amplitude ratio is not particularly sensitive to mode 
shape, permitting any reasonable fundamental type mode to be assumed. It 
must, of course, be ascertained, that the loading or forcing function is not 
of the type that would preferentially cause higher type mode-forms to be 
developed. With strain gages mounted at the center of the rod the second 
mode would not be picked up. Also, since the energy required to produce a 
mode with a given amplitude varies as the fourth power of the mode number 
it is apparent that a large energy input would be required to maintain vi- 
brations of any significant amplitude in the higher mode shapes. 

In order to present as clearly as possible the behavior of vibrating rods 
with elastic supports, the discussion has been limited to rods with equal 
elastic restraints at each end. The problem of vibrating rods with unequal 
end restraints is an interesting one, and will be presented at some future 
date. 
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SUMMARY 


This paper reports an experimental investigation of eccentrically loaded 
hinged-end steel columns bent about a minor axis parallel to the flanges. (1) 
In this respect it differs from, and supplements most existing experimental 
column investigations which are chiefly of two kinds: I-shaped columns tested 
(1) with bending applied about a minor axis parallel to the web, so that maxi- 
mum stresses appear only at the flange tips, or (2) with bending applied about 
a major axis parallel to the flanges, in which case flexural failure is normally 
accompanied by torsion. The present tests, in which bending is about a minor 
axis parallel to the flanges, preclude torsion while at the same time producing 
maximum stresses over the entire width of the flanges. 

Twenty-four eccentrically loaded and six concentrically loaded columns 
have been tested. In addition to initial yielding and ultimate loads, information 
has been obtained on load-deflection performance not only up to, but consider- 
ably beyond the ultimate loads, a matter which is of some consequence for 
plastic strength analysis of structures. No new theories or formulas for 
strength prediction are developed, but results are compared with predicted 
failure loads calculated by a variety of methods, including the secant formula, 
a plastic strength formula, and two interaction type formulas. 

In addition to these 30 new tests, the results of 28 tests on eccentrically 
loaded, elastically end restrained columns which have been reported previous- 
ly are compared in the Appendix with values calculated by the same two inter- 
action formulas. 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1792 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 

of the American Society of Civil Engineers , Vol. 84, No. EM 4, October, 1958. 
Asst. Prof., Dept. of Structural Eng., Cornell Univ., Ithaca, N. Y. 
Associate Prof., Dept. of Structural Eng., Cornell Univ., Ithaca, N. Y. 
Prof. and Head, Dept. of Structural Eng., Cornell Univ., Ithaca, N. Y. 
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INTRODUCTION 


Presently available experimental data on the strength of eccentrically- 
loaded two-flange columns, predominantly I-sections, are incomplete in two 
respects. First, there is considerable experimental data on the strength of 
two-flange sections eccentrically loaded to produce bending about an axis 
perpendicular to the flanges which, for the usual I-sections, is the minor axis. 
Jn this case maximum stresses occur at the tips of the flanges only, and ob- 
served ultimate loads appreciably exceed the strength computed by initial - 
yield formulas such as the secant formula. Second, in other tests bending was 
induced about the major axis of unbraced I-shaped columns. Such bending is 
almost always accompanied by twist, and test results cannot be interpreted by 
the usual formulas which exclude this effect. However, many structures in- 
corporate I-shaped columns bent about the major axis, but so braced against 
minor-axis buckling that twist does not occur. 

Therefore, available test data are neither wholly adequate for comparison 
with present column formulas, nor are they wholly applicable to the commonly- 
occurring condition in many structures where I-shaped columns are bent 
about the major axis parallel to the flanges and at the same time braced 
against buckling about the minor axis, effectively excluding twisting. Like- 
wise, twisting does not occur for box-shaped members bent about either axis. 

In order to overcome these deficiencies, hat-shaped sections, formed by 
connecting two rolled Z-sections flange-to-flange, were fabricated to produce 
a minor axis parallel to the flanges and were loaded eccentrically to produce 
bending about this axis. Consequently, the maximum compression stress was 
uniform over the entire width of the compression flange, and since bending 
was about the minor axis, there was no tendency of the entire column to twist. 
Therefore, results could be evaluated by methods applying to non-torsional 
flexural column action. 

In addition to failure loads, the behavior of such members beyond the ulti- 
mate load is of significance in relation to ultimate strength design, and for 
this purpose load-deflection observations were made far into the descending 
branch of the load-deflection curves. Width-to-thickness ratio of flanges was 
varied to investigate the effect of flange dimensions on this post-ultimate be- 
havior and on the ultimate strength. 


Specimens 


Test specimens consisted of a pair of Z-sections placed to form a hat- 
shaped cross-section, with outstanding flanges connected by batten plates and 
connected flanges stitch-welded intermittently. See Fig. 1. Thirty columns 
were fabricated and tested, using three different weights of rolled Z-sections, 
Three lengths of columns were made of each weight, and eccentricity ratios 
ec/r2 were varied from zero to 1.5. Specimen properties are given in Table 1. 


The following compression yield points were used in evaluating test results 
for the column specimens: 


1/4x 3 (3x 2 3/4, 6.7# 42,500 psi. 
1/4 x 4 (4x 3, 8.20 44,500 psi. 
1/2 x 3(3 x 23/4, 12.6 39,000 psi. 
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TABLE 1 


PROPERTIES OF SPECIMENS AND ECCENTRICITY RATIOS TESTED 
Z-Section L/r ec/r* b/t 


3x 2 3/4, 6.14 f 3.94 0; 0.25; 0.75; 1.50 10.8 
3x 2 3/4, 6.78 3.94 0; 0.25; 0.75; 1.50 10.8 
3x 2 3/4, 6.7# 3.94 0; 0.25; 0.75; 1.50 10.8 


4 x 3, 8.24% 1/4 4.82 36 0; 0.25; 0.75; 1.50 12.3 
4 x 3, 8.24% 1/4 4.82 107 66 0; 0.25; 0.75; 1.50 12.3 
k x 3, 8.24% 1/s 4.82 181 110.5 0; 0.25; 0.75; 1.22 32.3 


3x 2 3/4, 12.64 1/2 7.38 59 53 0.25; 1.50 5.4 
3x 2 3/4, 12.64 1/2 7.38 83 74 0.25; 1.50 5.4 
3x 2 3/4, 12.64 1/2 7.38 isi | (87 0.25; 1.50 5.4 


Notes: (1) A, r, ¢c, b, and t are based on nominal section properties. 
b = width, t = thickness of flange of Z-section. 
(2) Areas were spot-checked on a weight basis and nominal areas 


found to be sufficiently reliable for evaluation of tests. 


These were based on average values observed from compression stress- 
strain tests of short lengths of whole sections in the unannealed (as-received) 
condition, the results of which were compared for consistency with (a) simi- 
lar tests on annealed material to assess the influence of residual stresses, 
(see Fig. 2), (b) tensile coupon tests of material as-received, and (c) residual 


stress measurements following the method described in Lehigh tests, (3) Typi- 
cal residual stress measurements are given in Fig. 3. 


Testing Procedure 


The column testing arrangement, instrumentation, and centering procedure 
were similar to those described in previous papers except for end 
fixtures, (2,4) Briefly, the columns were supported on simple knife-edge 
fixtures parallel to the minor axis, were equipped with set screws to permit 
eccentricity adjustments, and were loaded so that maximum compression oc- 
curred in the outstanding flanges. 

Dial gages were located at mid-height and quarter-points of the column, 
and at the upper head of the testing machine. Tuckerman optical strain gages 
were placed at both quarter-points for centering, then at the mid-height only 
for testing. Four gages were used at each quarter-point, one pair on the 
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outer edges of the compression flanges and one pair on the edges of the 
tension face; at midheight, six gages were used: one pair each at inner and 
outer edges of the compression flanges and one pair at the outer edges of the 
tension face. 

For column centering, control was achieved with two aids: (a) computed 
charts of load-vs-midheight deflection for various eccentricities, based on 
elastic behavior and (b) theoretical yield loads computed according to the 
secant formula as a guide for limiting maximum centering loads and selecting 
load increments to be used in the tests. Final centering loads were about 60 
to 65 per cent of the computed yield loads. Centering procedure was as 
follows: (1) preset end eccentricities to desired value, (2) equalize lateral 
stress distribution by balancing strains on opposing quarter-point strain gages 
and shimming the column ends as required, (3) adjust end eccentricities 
slightly to produce midheight deflections as computed, and equal quarter-point 
deflections. It may be noted that deflection centering using dial gages is sig- 
nificantly more accurate than centering by use of strain gages. 

When centering was considered satisfactory, the strain gages were moved 
to midheight for the test. Strain readings were taken to the limit of range of 
the Tuckerman gages. Location and progression of yielding, as indicated by 
a white-wash stress coat, were recorded and particular attention was paid to 
flange buckling, noting when and if it did occur. Deflections were read 
throughout the test, straining being continued beyond the ultimate load until 
the equilibrium load had dropped to about 65-75 per cent of the ultimate load. 


Behavior Under Load 


Deflections agreed satisfactorily with those predicted by elastic analysis 
nearly to the computed yield load. First yielding occurred at loads somewhat 
less than the predicted values, attributable to the high residual stresses ob- 
served. Ultimate equilibrium loads were usually slightly above the yield loads 
predicted by secant formula. (See Table 2). 

Behavior beyond the ultimate load followed a consistent pattern. Loads 
dropped off quite slowly for columns with large eccentricity as opposed to 
more rapid loss of capacity for concentric columns and those with small 
eccentricity. Likewise, loads dropped off more quickly for low slenderness 
ratios whereas for the most slender columns loads held at or near the maxi- 
mum over a considerable range of deflections (see Fig. 4). At ultimate loads 
midheight deflections amounted to 0.2% to 1.15% of the length, which is a 
moderate magnitude, considering that deflections at design loads would be a 
fraction of these values. 

Column strength was not affected by local flange buckling, since when it oc- 
curred at all, it did so well beyond the ultimate load. Likewise, no influence 
could be detected of flange buckling on the post-ultimate load-deflection be- 
havior. It had been conjectured that for large b/t-ratios, early flange buckling 
at loads just past the ultimate might lead to quicker collapse, i.e. a steeper 
load-deflection curve, than for stockier sections, but no such trend of signifi- 
cant magnitude was observed within the range of dimensions investigated, i.e. 
within w/t from 4.4 to 11.3, where t = thickness, w = projection of outstanding 
flange. 

The point of initial yield, not reliably indicated by the white-wash stress 
coat, was determined by a semi-empirical method described in a previous 
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paper. (2) This consisted of computing the midheight bending moment from 
the load and the corresponding total midheight eccentricity, computing there- 
from the maximum fiber stress, and plotting this stress against the observed 
strain to obtain a linearized stress-strain curve. Initial yield load occurs 
where this plot first deviates from linearity. Initial yielding so determined is 
indicated by the dashed curves of Fig. 4 and is seen to occur typically at loads 
but slightly smaller than the ultimates. 


Ultimate Loads: Test vs. Prediction 


Experimental ultimate loads have been compared with predicted values by 
the secant formula, by the Bijlaard method, (2) and by two interaction formu- 
las. For hinged-end columns the ultimate load by the Bijlaard method is 
quite simply obtained since the effective length is known. A shape factor 
HM = 1.4 was used as suggested in the previous publication (2) for bending about 


an axis parallel to the flanges. The mean eccentric buckling stress is ex- 
pressed by the formula: 


(t)e is the slenderness ratio for which the buckling stress is equal to the 
proportional limit (assumed to be 10,000 psi less than the yield point). 


(+) - 3x 10° 
10,000 
and where the reduction factor % which accounts for eccentricity is obtained 
as follows, where () = 157 - 0.0014 o,., 
Tr y 


for 200 > 


(¢) 


1+ (21000 - 0.23 ¢,) — 


L L 

where = 50,000 7% + (5) 

for - 30 > > 20 
rh 
e K. a 

Ke 370 - 0.005 


ASCE HINGED STEEL COLUMNS 1792-9 


L L 
for (F) > (=) - 30 interpolate linearly between ¥ values 


for - 30 and 


cross-sectional area of column 
eccentricity of load 


core radius 


= 
length of column 
¥ r = radius of gyration of cross-section = 
= division factor for eccentricity 
H = shape factor q 
0, = mean equilibrium stress 2 
cES = Engesser-Shanley buckling stress (by Bijlaard method) 
0 ys = yield point stress 7 


Then, the predicted ultimate load is 


Put = xA 


Table 2 compares measured ultimate loads with those calculated by the 
secant formula and the Bijlaard method. The latter is seen to be slightly more 
conservative than the former, but both are satisfactorily accurate in predict- 
ing the capacities of these columns. 

From Table 2, the Bijlaard method is seen to give good correlation for 
columns of short and intermediate lengths, while the predicted ultimate loads 
for the long specimens are consistently below the experimental values. Thus 
it would appear that the value for the effective proportional limit which has 
been assumed in determining the Engesser-Shanley buckling stress is too 
conservative when applied to these double Z-section columns, 

From the same table it can be seen that the secant formula gives very 
satisfactory correlation with experimental values for ultimate loads for all 


slenderness-ratios. 
Failure loads have also been evaluated according to two types of interaction % 
formulas. The simpler type is the one current in several widely used design fe 
codes, namely “Interaction Formula Type I”: 
P M 7 


where P = actual load at failure under combined action 


M = bending moment at failure under combined action 


P' = axial load at failure under concentric compression, the 


Engesser-Shanley load (when determined by the Bijlaard 
method, P' = d,.pgA) 
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ec Pye Ratio 
Test Secant 

kips kips 1) /(2 

4g 0.25 «122 123 0.992 
69 0.25 110 ll 0.991 
108 (0.25 Th 1.007 
4g 0.75 89.3 85.5 1.044 
69 0.75 77-6 76.8 1.011 
108 75 57.2 56.7 1.008 
4g 81.50 65.9 60.7 1.085 
69 1.50 58.4 54.8 1.066 
108 1.50 42.5 1.016 
36 0.25) 164 1.048 
66 0.25 143.8 143 1.006 
110.5 0.25 87.8 90.1 0.974 
36 «600.75 «123.2 15 1.070 
66 0.75 100.1 99.0 1.012 
110.5 0.75 66.2 68.5 0.967 
36 4861.50 8h.2 81.0 1.039 
66 1.50 71.0 71.2 0.997 
110.5 1.215 58.1 57.9 1.004 
53 0.25 211 1.015 
T 0.25 188.6 187.5 1.005 
122 1.002 
53 1.50 117.2 103 1.138 
Tm 1.50 101.2 93.0 1.088 
117 1.50 76.1 71.0 1.071 
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TABLE 2 
ULTIMATE LOADS - TEST VS. THEORY , 
Specinen Ratio 
in. kips) (2 3 
x 3 118 1.034 
1/4 x 3 105 1.047 
x3 66.4 1.129 
3 87.8 1.017 
x3 76.5 1.01% 
x 3 1.158 
1/4 x3 64.5 1.022 
1/4 x 3 56.2 1.039 
1/h x 3 36.1 1.134 
xb 162 1.060 
x 135 1.065 
1/4 81.2 
exh 121.5 
xh 9.0 
1/4 xh 61.7 1.073 
xb 89.7 0.939 
xh 72.7 0.977 
xh 52.4 1.108 
1/2 x 3 200 1.071 ‘i 
1/2 x 3 173 1.090 
1/2 x 3 109.5 1.117 y 
1/2 x 3 109.5 1.070 
1/2 x3 91.8 1.103 
1/2 x 3 64.4 1.182 ® 
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M' = bending moment at failure under bending only. For simple 
bending, it is the plastic hinge moment, or if lateral buckling 
occurs, the critical lateral buckling moment. 


The “Interaction Formula Type II”, is: 


M 


+ M (1 - P/Pp) P) Pr 
where Py =the Euler elastic buckling load. 


Comparisons are given in Table 3 and Figs. 5 and 6. 

Inspection of Fig. 5 reveals that the “simple”, widely used interaction 
formula of Type I is unconservative. Of all eccentric tests, only one reached 
an ultimate load barely exceeding that predicted by this formula. In all other 
cases the formula overestimates the actual strength, frequently by 15 to 20 
per cent and more, and also results in rather wide scattering. On the other 
hand, Fig. 6 shows that the interaction formula of Type I, is safe in nearly all 
cases, is never excessively conservative, and results in rather narrow scatter- 
ing. 

This comparison is somewhat distorted because the values for P' (for zero 
eccentricity), computed by the Bijlaard method from the Engesser-Shanley 
stress O.~g (see above), are somewhat too conservative in that they over- 
estimate, for these sections, the influence of residual stresses in lowering the 
effective proportional limit. This was mentioned previously, is apparent from 
Table 4, and is reflected in the interaction plots where for M = 0 all test re- 
sults fall above 1 (one). If a less conservative approach had been used in com- 
puting P', interaction formula I would have been more unsafe, while formula I 
would have resulted in a straight line very close to the middle of the narrow 
scattering band. 

A statistical study of the correlation of test data with the various theories 
is summarized in Table 5. For the secant formula and the Bijlaard method, 
the ratio of the experimental ultimate load to the predicted value was used as 
the “observation”; the sum P/P' + M/M' was used for interaction formula I, 


P M 
and the sum for interaction formula II. Arithmetic means, 


standard deviations, and coefficients of variation have been computed for an 
assumed normal distribution, since the number of tests was too small to justi- 
fy further refinement. 

The results of the statistical study, as given in Table 5, show quite definite- 
ly that for these columns, bent about an axis parallel to the flanges, and with- 
out tendency to torsional instability, the secant formula gives the most satis- 
factory correlation of predicted and actual ultimate loads. Next in agreement 
are the interaction formula of Type II and the Bijlaard method, both of which 
show satisfactory and practically identical coefficients of variation, and both 
of which err slightly and systematically on the conservative side, for reasons 
previously discussed. The customary interaction formula, Type I, shows the 
largest coefficient of variation, i.e., the largest scattering range, is the only 
one which consistently errs on the unsafe side, and must be regarded as 
unsatisfactory. 
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TABLE 3 
INTERACTION FORMULAS 


I 
P 
P' 


Specimen L/, 


1/4 x 3 
1/4 x 3 
x 3 
1/4 x 3 
1/4 x 3 
1/4 x 3 
1/4 x 3 
1/4 x 3 
1/h x 3 
xh 
1/4 x 4 
exh 
exh 
1/4 
1/2 x3 
1/2 x3 
1/2 x 3 
1/2 x 3 
1/2 x 3 


1/2 
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tx bd (1) + (2) (1) + (3) 4 

49 0.25 0.802 0.156 0.958 0.210 1.012 ‘ ™ 

69 0.25 0.769 0O.1b2 0.910 0.257 1.026 7 

108 0.25 0.768 0.096 0.864 0.412 1.180 

49 0.75 0.587 0.343 0.930 0.423 1.010 

69 0.75 0.542 0.298 0.840 0.438 0.980 
108 0.75 0.586 0.220 0.806 0.528 

491.50 0.433 0.508 0.9h1 0.589 1.022 

69 1.50 0.408 O.4kg 0.857 0.591 0.999 x 

108 1.50 0.443 0.333 0.776 0.598 1.041 j 

36 0.25 0.850 0.17% 1.024 0.208 1.058 

66 0.25 0.778 0.145 0.923 0.263 1.041 

110.5 0.25 0.763 0.089 0.852 0.376 1.139 _ 

36 0.75 0.610 0.375 0.985 1.034 7 

66 0.75 0.541 0.305 0.846 0.983 

110.5 0.75 0.576 0.202 0.778 1.052 

36 1.50 0.417 0.514 0.931 0.557 0.97% 

66 1.50 0.364 0.432 0.816 0.555 0.939 be 

110.5 1.215 0.505 0.287 0.792 0.579 1.084 Me 

53 0.25 0.828 0.1h9 0.977 0.206 1.034 2 

TH 0.25 0.783 0.132 0.914 0.248 1.031 

U7 0.25 0.779 0.085 0.864 0.385 1.164 ; : 

53 1.50 0.453 0.488 0.941 0.577 1.030 Y 

Te 1.50 0.420 0.422 0.842 0.566 0.988 

1171.50) 0.485 (0.317 0.802 0.615 1.100 
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TABLE 


CONCENTRIC SPECIMENS 


(1) (2) (3) 4 (5 
Specimen L/. Equilibrium  Meximm# Theor. (1)/(3) (2)/(3) 
Kips P, Kips P', Kipe 


1/4 x 3 156.8 162.6 152.2 1.030 1.068 
ifs x3 163.4 166.5 143.0 
1jex3 10 102.0 104.1 98.0 1.041 1.062 


ijexh 3% 203.5 208.5 202.0 1.007 1.032 
ijexh 6 189.0 193.7 185.0 1.022 1.047 
iaexh 110.5 119.7 121.2 115.0 1.040 1.054 


* Highest Stable Load 


* Maximm load reached, but load not stable 


Note: Stiffening effect of batten plates and end blocks 
has been neglected. 


TABLE 5 
COMPARISON OF VARIOUS FAILURE THEORIES 


Type of Range of Arithmetic Standard Coefficient 
Observation Observations Mean Deviation of Variation 


Secant Forma 0.967 to 1.138 1.027 0.0398 3.87% 
Bijlaard Method 0.939 to 1.182 1.065 0.0562 5 .20% 
Interaction Type I 0.776 to 1.02h 0.882 0.0658 1.47% 


Interaction Type II 
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Discussion 


It is seen that the ultimate loads of these columns were predicted with al- 
most equal accuracy by the secant formula (i.e. based on elastic calculations), 
by the Bijlaard method (based on plastic calculations), and by the interaction 
formula of Type II. The only method which proved definitely unsatisfactory is 
the interaction formula of Type I. 

In regard to the three satisfactory methods, the following may be said in 
regard to their validity for columns of types other than those tested: 

The secant formula (or any determination based on an initial-yield criteri- 
on) is accurate only for bending about an axis parallel to the flanges, and pro- 
vided that such bending is not accompanied by torsion. It has been shown in 
Ref, 2 that in the case of bending about an axis parallel to the web the initial- 
yield criterion underestimates the column capacity considerably, up to 100% 
in those tests. On the other hand it is well known (see e.g. Ref. 5) that for 
columns bent about a major axis parallel to the flanges, torsion can reduce the 
ultimate load considerably below that predicted by an initial-yield criterion. 

As to the Bijlaard method, the combined results of the present tests and 
those reported in Ref. 2 indicate that this method is satisfactorily accurate no 
matter whether bending is about an axis parallel or perpendicular to the 
flanges, and no matter whether the columns are hinged or elastically re- 
strained. The principal limitation is that it does not apply in the case of 
flexural-torsional buckling. That is, in case of application to bending about a 
major axis parallel to the flanges, it is restricted to the case where torsion 
is prevented as by appropriate minor-axis bracing. 

The interaction formula of Type II has been proposed by various investi- 
gators, (6,7) In the present paper and its Appendix it is shown that the results 
of this formula are satisfactorily accurate for bending about axes either paral- 
lel or perpendicular to the flanges when torsion is prevented, and for end- 
restrained as well as for hinged columns, In addition, however, it has been 
shown previously by others(7,8) that this formula is also satisfactory for 
aluminum columns subject to flexural-torsional as well as purely flexural 
failure. Finally, by means of an extensive series of tests, C. Massonnet and 
F. Campus(9 have shown that for steel columns also this formula is likewise 
satisfactory both for flexural and for torsional-flexural failure, including the 
case of unequal as well as opposite end moments. In the case of torsional- 
flexural failure (major-axis bending of unbraced I-shaped columns) the criti- 
cal moment for lateral beam buckling is substituted for M', and for unequal 
end moments an equivalent moment is used for M which is obtained from a 
very simple formula. 

To summarize: While it is possible to devise satisfactory formulas for 
various special cases of eccentrically loaded columns, of those so far pro- 
posed the only simple equation which seems to apply satisfactorily to all 
situations is the interaction formula of Type I. 


CONCLUSIONS 


1. Although flange-projection-to-thickness ratios ranged from 4.4 to 11.3, 
flange buckling occurred only in the plastic region, did not affect column 
strength, and had no discernible effect on the behavior of the columns 
past the ultimate load, in that region of large deflections which is 
important in connection with ultimate strength design. 
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2. For these columns bent about an axis parallel to the flanges the secant 
formula showed best correlation with experimental data. It is conclud- 
ed that for columns bent so that the entire flange width is stressed uni- 
formly, the ultimate load is practically identical with that for which the 
elastically computed flange stress becomes equal to the yield point, 
provided torsional-flexural failure is precluded. 

3. The Bijlaard method also was satisfactory for these particular Z- 
section columns, Correlation was not quite as good as with the secant 
formula, but improvement would be possible by adjustment of the as- 

sumed proportional limit stress, in order to reflect more nearly the 

small influence of residual stresses in these tests. 


The simple interaction formula, = + x = 1 is definitely unconservative 


for eccentrically loaded columns. For all but one of 24 specimens this 
formula predicted larger strengths than were obtained by test. 


5. The interaction formula, (i - P/Pp) 1 showed very satis 


factory correlation with test data. Scattering was small and only six 
specimens failed at loads lower than the predicted values, the largest 
error in this direction being 6 per cent and the rest of the “unsafe” 
observations erring by 2-1/2 per cent or less. 

6. Evidence presented in the Appendix shows that the better interaction 
formula is also satisfactory for end-restrained columns, although it 
does not provide a significant simplification in this case. 

7. Reviewing the evidence of other investigators it is concluded that only 
the interaction formula of 5, above is equally satisfactory for all in- 
vestigated cases of eccentric columns: those of steel or aluminum, 
hinged or restrained, with equal or unequal eccentricities, with bending 


parallel or perpendicular to the flange, and failing flexurally or by com- 
bined torsion and flexure. 
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Appendix 


End-Restrained Columns en 
An extensive experimental and theoretical investigation on elastically end- sa 
restrained, eccentrically loaded columns was reported in Ref. 2. Column 
specimens consisted of square-bar sections and of I-sections bent about the iy 
minor axis, parallel to the web. It was shown that the Bijlaard method, as ap- a r 


plied to these restrained columns, predicted the failure loads very accurately. 
On the other hand, it was shown that the measured ultimate loads exceeded 
greatly (up to 100%) the calculated incipient-yield loads of these columns (i.e. 
the loads calculated for these restrained columns on the same basis on which 
secant-formula loads are calculated for hinged columns.) 

In view of the simpli¢ity of interaction formulas, and in view of the versa- 
tility of the interaction formula of Type II as explained in the body of the 
paper, it seemed of interest to evaluate the results of that investigation in 
terms of these two interaction formulas. This is done in the present appendix. 

As indicated in Ref. 2, elastically-restrained eccentric columns of length 
L and eccentricity e are successfully analyzed by replacing them by equivalent 
eccentric hinged-hinged columns of length cL and with effective eccentricity 
€e. The correct determination of cL and eg involves the calculation of the 
plastic reduction factor 7, which is computed by successive approximation. 

Figs. 7 and 8 were constructed on this basis. That is, in these figures 


P is the ultimate test load 


ad is the concentric Engesser-Shanley buckling load determined from 
Eq. (20) of Ref. 2 for the equivalent column of length cL, where this 
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length is calculated on the basis of the appropriate plastic reduction 
factor 7. 


is the portion of the total bending moment Pe which at failure, is re- 
sisted by the column, i.e. Mp = Peg, where the equivalent eccentricity 


e, is again determined by means of the plastic reduction factor n, 
as indicated in Step 5, p. 1086 of Ref. 2. 


is the plastic hinge moment of the member about the axis about which 
bending was applied (the minor axis for the I-sections), computed by 
the simple plastic theory from the dimensions of the section and the 
yieid point obtained from the compression tests. 


is the elastic Euler load of a column of length c'L, where c'L is de- 
termined for fully elastic conditions, i.e. for 7 = 1. (Methods for 
this determination are given in Ref. 10). 


The above method is intrinsically consistent with the analytical approach 
developed in Ref. 2. Its practical disadvantage, in connection with an inter- 
action formula, consists in the fact that use is made of the plastic reduction 
factor 7, which must be determined by successive approximation. Hence, the 
use of an interaction formula, if successful, merely replaces Step 7, p. 1088 
of Ref. 2, which happens to be the simplest and most direct of the various 
steps there given. For this reason interaction formulas have little real ad- 
vantage in this case. Their validity is explored in this Appendix chiefly to 
verify whether they apply to restrained columns. 

Fig. 8 indicates that the use of the more accurate interaction formula, Eq. 
Il, results in a grouping of test values close to and above the straight line 
representing this equation. Hence, if this equation is used in connection with 
the correct plastic reduction factor 7, a safe and reasonably accurate de- 
termination of failure loads results. 

Fig. 7 indicates that the simpler interaction formula, Type I, used in con- 
junction with the correct plastic reduction factor n, gives significantly 
unconservative determinations for a sizeable number of the tests. This is in 
agreement with the findings of other investigators and with the evidence pre- 
sented in the body of the present paper. 


+2 
ASCE 
3 
7 
M 
4 
* 
P 4 
| 
< 
q 
ve 
A 
~ 
>. ‘ha 
q 
ae 
4 
in 
a 


= 
be 
a 


Paper 1793 October, 1958 EM4 


Proceedings of the American Society of Civil Engineers 


Journal of the 


ENGINEERING MECHANICS DIVISION 


APPROXIMATE BUCKLING LOADS OF OPEN COLUMNS 
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SYNOPSIS 


A general equation for buckling loads of open thin-walled columns is pre- 
sented in this paper. The equation is applicable to any boundary conditions. 
The coupling effect between torsional buckling and bending bucklings is con- 
sidered. 
as two double coupling cases without introducing significant errors. Equations 
are also given for columns restrained by flexible skin. 


It is found that a most general case of triple coupling may be treated 


NOTATIONS 


Cross-sectional area 


Uniform torsion constant (corresponding to the Saint-Venant Theory 
of Torsion) 


Warping constant, = he (D - w,)? dA 
Average warping of a section due to unit twist 
Modulus of elasticity 

Modulus of elasticity in shear 

Moment of inertia 

Polar moment of inertia 

External end moment 

Resultant end thrust 


Strain energy 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1793 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 4, October, 1958. 


1. Prof. of Civ. Eng., Imperial College of Engineering, Addis Ababa, Ethiopia. 
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W Work 


» 4 Dimensionless function of x, representing mode shape 

Q Length of column 

m Mode shape factor 

Ws Warping at any point of a section due to unit twist 

Vv, Ww Translations in the y, z directions respectively (referred to the 


shear center without subscription 


x, y, Z Coordinate axes through shear center, or the corresponding co- 
ordinate values 


Y Coupling factor 

8, 2 Centroidal coordinate axes (7 and ¢ are parallel to y and z re- 
spectively) 

K A section property 

to) Angle of twist 

Unit compressive stress 


The following symbols are used for subscripts 


Cc Centroid 

M, N Line oi attachment to skin 

Shear center 
Q Point of application of resultant thrust c 
m Modified 


INTRODUCTION 


A prismatic member of thin-walled open section may have small rigidity 
against torsion. When used as a column, it may buckle in torsion, or in 
coupled bending and torsion. In general, the solution of the critical load, at 
which the buckling occurs, involves three simultaneous differential equations, 
and usually is difficult to solve except in the case of a column with simple ¥ 
supports at both ends. Further difficulty arises when a mixed boundary is in 4 
question, such as an end condition providing a simple support for torsion but + 
a clamped support for bending. Using the Rayleigh-Ritz method these compli- 4 


cations can be avoided, and the solution can be extended to cover all boundary 
conditions. 


Energy Consideration 


The Rayleigh-Ritz method is based on the principle of conservation of 
energy of a system. Let us consider a column of arbitrary thin-walled open 
shape as shown in Fig. 1. C is the centroid, and 0 is the shear center of the 
section. 7 and ¢ are the centroidal principal axes corresponding to the 


= 
359 


ASCE OPEN COLUMNS 1793-3 


2 
7 
a 
a 
: 
all 
Q 
SSN 
fi 
Fj 
ig. | 
M 
4 q 
oO 
8 
‘ 
Q YC 
a / 
4 
Pp 
My= P (29 - 
1 
j 


1793-4 EM 4 October, 1958 


minimum and maximum moments of inertia. We take the coordinate system 
with the origin at the shear center 0 of an end section, and with the y and z 
axes parallel to 7 and £ respectively. In this way, the x axis will coincide 
with the locus of shear center when the column is undeformed. 

The column is loaded at each end with a thrust P through point Q. For 
convenience of investigation, P will be replaced by an equivalent force system 
as shown in Fig. 2, where P now passes through C, and Mn, M;, are moments 
about the and _ axes, respectively. 

Any displaced position of a cross section can be completely defined by the 
translation components v and w of the shear center of the section (the dis- 
placed shear center is indicated as 0' in Fig. 3) and a rotation ¢ around 0'. 
Owing to the end moments Mn, M+, there are some initial deformations be- 
fore the buckling. These initial deformations will be disregarded in the 
following study of energy transformation associated with the buckling only. 
More precisely, the components of displacement, v, w and @ in the discussions 
to follow are the additional displacements measured from the initially de- 
formed shape just before the instant of buckling. 

The additional strain energy due to additional displacements associated 
with the buckling may be expressed as follows:2 


where the primes represent the order of derivation with respect to x. It is 
seen that Eq. (1) does not include the energy contributed by the direct shear 
(distinguished from torsional shear) and the additional shear due to warping. 
This is permissible provided that the column is not extremely short.3 

The increment of energy given in Eq. (1) results from the additional work 
done by the end thrust during buckling. The end thrust on an elemental fiber 
has an intensity of 


M M 
)- 


(1) 


A 


The sign convention used in Eq. (2) is consistent with Figs. 1 and 2, anda 
positive sign corresponds to compressive stress. When buckling occurs, the 
two ends of the fiber move toward each other by an amount equal to 


(2) 


2. Ref. (1), Eq. (29), after setting the product of inertia equal to zero. Refer- 
ence numbers refer to similarly numbered references in bibliography at 
the end of the paper. 

3. For avery short column, failure would usually be a local crippling type 
with original shape of cross-section distorted. See Ref. (2). 
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(4) 
Observing that 


A 
A 
A 
A 


A 


and denoting 


k -{ ( y= aa =f y ( y2#2*) da = dy 
4 7A A 


4. By definition, D is the average value of wg of a section, and so the relation 


is self evident. 
5. See Ref. (1), Appendix III. 


z dA = 2, A 
| y dA = Yo A 4 
( y?+2°) di=dy 
(y- Yo ) dA = 0 
( 2-2.) dA #0 q 
A A 1 
an 
yliy-y, )da=fly-y di = 
(D-w,) 
A 
| (y-yo )( D-w, ) di 
A A 
= 0 5 : 
dt - 2, Jy 
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we obtain 


Mn My (2 
+2 w"dx ~2 \v vidx| 
P Jo 
Substituting 
Mr = P (zg - zc) 
M; = -Pp (yg yc) 


Eq. (4') may be written as oa 


2z +24, +f (ow)? dx 


where 
Ky k 


In general, we may express 


U 


‘Vv 
Wo 


where vo, wo, $09 are amplitudes, and Xy, Xy, X¢, are dimensionless functions 
of x, representing the buckled shapes of the first mode® in three directions 
respectively. Substituting (5) into Eqs. (1) and (4''), we get 


AU=: | vo EI 
GC aj 4 EC (6) 


6. The lowest mode order corresponds to the smallest of the eigen values of 


the critical load which is the only one of practical importance. .. 
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AU= 53K (6) 


{ 2 24m 4\2 


4+2v, 4 x| (7) 


= Sh (7) 


where & and § represent the expressions inside the brackets of (6) and (7) re- 
spectively. 


Critical Loads for Coupled Buckling 


The critical load at which buckling will occur may be obtained by equating 
(6) and (7), which gives 


(8) 


It is seen that both @ and # contain arbitrary amplitudes Vo, Wo, oo. The 
true values of these amplitudes should render a minimum value for P. This 
condition is satisfied when 


Applying (9) to (8) and multiplying the resultant expression by 8 we obtain 


1793-8 

and 

or 

0 Yo 

oP 

—— = 0 (9) a 
2 Wo 

x 


Substitution of the expressions for o@ and B into Eqs. (10) leads to three simul- 


taneous equations for vg, wo, and $9, namely 
+2(4,- ydx]=0 --------(a) 
+2(_- = -------. (b) 
After rearrangement of terms, we get 
[E Yax -P [P (x, x,42| 5, 
=P(4- 4 (a’) 
=P(Ze dx --------- (b) 
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The values of vg, Wo, and do become indefinite (mathematical interpretation 
of instability or buckling) when the following determinate vanishes 


9 2 


0 
PI 


Q 
(x 

-P 


(11) 
= 0 


Eq. (11) is a cubic equation for p. The terms containing Cx Xy AX and 

j z ‘dx demonstrate the interaction (or coupling) between bending and 

torsio bucklings. The equation gives three positive roots corresponding to 

three possible ways of coupling. However, the smallest root is the only one 

of practical significance. 


Let us now consider a special case, in which yo = ZQ = 0 (resultant thrust 
P passes through the shear center 0). Steseninding We, (11), the interacting 
terms will no longer exist, and the three roots are evidently 


Ely 
( 
2 


The first two roots in (12) are critical loads corresponding to the noncoupled 
bending (or Euler) type, and the third root is the critical load corresponding 
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to the noncoupled torsional (or Wagner) type of buckling. We shall call Py, 
Pw and Pg the noncoupled critical loads in the v, w and @ directions re- 
spectively. 

Generally speaking, the coupling effect will alter the mode shapes in some 
degree, if the noncoupled bending mode shape are different from the non- 
coupled torsional mode shape for the specific boundary conditions. However 
such an alteration is usually small, if there is any, because the boundary con- 
ditions remain the same. In most cases, the noncoupled bending modes and 
the noncoupled torsional mode are identical,? then there will be no alteration 
of either shape. Even if there is some alteration of shapes, Eqs. (12) may 
still be considered as good approximate relations between the coupled mode 
shapes and the noncoupled critical loads. 


Using relations (12) in Eq. (11), and after expanding the determinant and 
simplifying, we obtain 


+P (Py Rp Fy Fa Fp 


(13) 


where 


2 
= a mod hape factor 


, 
= a moce shape fa 


Playa 


(14) 


When the noncoupled mode shapes in bending and torsion are identical (i.e. 
Xy = Xw = Xo), it is seen that myg = mwg = 1. And, 


+P(P, Fe + Fu * (15) 


is an exact equation for P for this case. But when X,, Xy and X¢ are not 


identical, myg and my¢ are less than unity (in view of Schwarz’s theorem of 
inequality). Their values can be approximated by substituting the noncoupled 


7. See Tables I and II in the following paragraph. 
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mode shapes (use Tables I and II in the following paragraph) in relations (14). 
In that case, Eq. (13) is to be used and the result is an approximate one. 
Fortunately, a slight change of the values of m affects very little the smallest 
root of Eq. (13) which is the only root important to engineers; therefore the 
approximate formula (13) is usually accurate enough for practical purposes. 

Eq. (15) has been given by several previous investigators in different no- 
tations and forms.8 They were derived for columns of simple supports at 
both ends with centroidal thrust. Here, we have proved that the same equation 
(or similar equation if Eq. (13) is used) may be applied to a column of other 
end conditions with a resultant end thrust at any arbitrary point. The only 
precautions are: (1) Py, Py, Pg represent the noncoupled critical loads 
corresponding to their respective boundary conditions; (2) A modified polar 
moment of inertia J, is used instead of the polar moment of inertia with re- 
spect to the shear center. 

Some conclusions drawn by the previous investigators for columns of 
simple supports and centroidal thrust9 may also be applied as well here. The 
smallest root of Eq. (15) or Eq. (13) is smaller than any of the noncoupled 
values Py, Py and Py, therefore disregarding the coupling effect is uncon- 
servative. However, if one of the noncoupled P’s is much smaller than the 
others, we can practically take that value as the smallest root, and the solution 
of the cubic equation (13) or (15) becomes unnecessary. For example, if 
Pg << Py, < Py, we shall find that the smallest P=Pg. The buckling is said 
to be torsionally predominant. On the other hand, if Py < Py << Pg, we may 
take the smallest P= P,. The buckling is said to be bending predominant 
(more precisely, in the v direction). If the three noncoupled P’s are very 
close to each other, the coupling effect will be large. 


Noncoupled Critical Loads and Mode Shapes 


For application of Eqs. (13) to (15) knowledge of the noncoupled critical 
loads and sometimes the noncoupled mode shapes is necessary. The non- 
coupled critical loads in bending (usually called Euler’s loads) are commonly 


known. They are given with the mode shapes in Table I corresponding to their 
boundary conditions. 


TABLE I 


Noncoupled Buckling In Bending 


Boundary Conditions Mode Shape (Xy or X,,) 


Tx 

1 - 


4.4.93 x x 
Simple - Fixed sin——>—_ - = sin 14.493 


Q g 


Tx 
Fixed Free 1 - cos 
(Cantilever) 2% 


Simple = Sinple sin 


Fixed = Fixed 


“& 


8. See, for examples Ref. (1), Eq. (20); Ref. 3, Eq. (84) and Ref. 4, Eq. (61). 
9. See, for example, Ref. (4). 
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The noncoupled torsional buckling was first investigated by H. Wagner, 10 
The differential equation for this type of buckling is 


x PJ, ” 
ECF ( GCu) $ = 0 (16)11 
&§ Its general solution is found to be 


t+ + Cy 


cA EC, 


There are four boundary conditions at two ends of a column, from which we 
can determine the valué A and three arbitrary constants in (17). One of the 
constants remains arbitrary. 

A specific interpretation for the boundary conditions in torsion is needed. 
A simply supported end means the column is prevented from twisting about 
any longitudinal axis, but is free to warp (or warping stresses are zero). 
When written in mathematical language, the two conditions correspond to 


By a fixed end in torsion, it is implied that both rotation and warping are pro- 
hibited, and 


are the mathematical expressions. For a free end, there is no restraint of 
any kind, and the following conditions must be satisfied 


= 0 
- 0 


The second condition means that the total torque about any longitudinal axis at 
the end should vanish. 12 

Applying the above-mentioned boundary conditions to (17), we obtain, simi- 
lar to Table I, the noncoupled critical loads and the noncoupled mode shapes in 
torsion as shown in Table II. 

In Table II, the interesting result is observed, that for a cantilever column 
the mode shape in torsion is zero, yet there is a definite value for the critical 
load. The fact is the free end offers no means to counteract any twisting, and 
the instability is probably of a local nature. 


10. Ref. (5). 
11. As stated before, noncoupled torsional buckling can only occur when the 
resultant thrust passes through the shear center, and Jp, = Jog. However, 


for convenience of later substitution, the notation J», is purposely retained. 
12. See Ref. (5), Eq. (5a). 
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TABLE II 


Noncoupled Buckling In Torsion 


Boundary Condition Mode Shape ( Xp) Pe 


2 
Simple -Simple GCuA. ECWA 
a 
im Fi i ; GCuA T CWA 


Fixed-Free 
(Cantilever) 


Some Comments on Eas. (13) and (15) 


It is noted from Table I and II that the noncoupled bending modes and 
torsional modes for the same boundary conditions are identical except in the 
case of a cantilever. As a result, Eq. (15) may be used most of the time, ex- 
cept in some mixed combinations such as an end condition providing a simple 
pgp for torsion but a clamped support for bending.13 For those cases, 


ase and Myq@ will have to be actually introduced and the approximate equation 
(13) used. Ad the resultant thrust P passes through the z axis, i.e. yg = 0, and 
consequently ”y¢ = 0, Eq. (13) is then reduced to 


It is evident that Eq. (18) has an independent root P = P,, which is the non- 
coupled critical load in bending in the w direction. The remaining two other 
directions of buckling, v and@ are coupled together. The coupled critical 
loads are represented by the roots of the following equation 


P(P,+ Fp )+ Fe (19) 


This special case is sometimes called the “Double Coupling” compared with 
the most general case of “Triple Coupling” when all three directions of buck- 
ling are coupled together. If the resultant thrust P passes through the shear 
center 0, i.e. yg = 0, zg = 0, Eq. (13) is reduced to 


(P- Py) Pu Fe) =0 (20) 
The roots are already given in (12) and they are all independent of one another. 


13. This is possible because the interpretation of boundary conditions in 
torsion is completely independent of that of bending. 
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Approximate Solution for a Triple Coupling by Two Double Couplings 


A study of Eq. (13) reveals that Yw@ has very little effect on the v di- 
rection. As shown in the last paragraph, Eq. (13) reduces to Eq. (18) when 
Ywo = 0. Similarly, when = 0, Eq. (13) reduces to 


(P-P.)| P(R,+ Py) > | =O 


A close approximation for the critical load P may be obtained by disregarding 
first wo and then Yo» one at a time, i.e. by solving only the quadratic 
portions of Eqs. (18) and (18'), The smallest of all the four roots may be 
taken as the critical load P sought. This procedure avoids the trial and error 
solution of a cubic equation. 

It should be noted, however, that the above-mentioned approximation is an 
upper-bound one. Fortunately, the error associated with the approximation 
rarely exceeds 2% for a column of ordinary proportions. In the most extreme 
cases, the error may reach about 10%. At any rate, the approximate value 
will provide a starting point for the trial and error solution of a cubic 


equation, if a very accurate result is desired. Re 
Column Attached to a Thin Sheet " 
The open section columns are often used as stiffeners in aircraft structures et 
and attached to the skin. The skin provides very large resistance on the line ba 
of attachment against movement in the plane of skin, but practically no re- - 
sistance perpendicular to the skin.14 We can idealize the problem as such = 
that the movement parallel to the plane of skin is zero on the line of at- = 
tachment. 
As shown in Fig. 4, the origin of the coordinate axes is still taken at the e 
shear center of an end section, but the y- and z-axes are now taken as perpen- Fe. 


dicular and parallel to the flexible skin. The 7 and ¢ are still centroidal 
axes and parallel to the y and z respectively, but are not necessary in the 
principal directions. Q is the point of application of the resultant thrust and 
N indicates the position of the line of attachment. 


Referring to these coordinate axes, the increment of strain energy and the 
additional work done by the end thrust are, respectively 


+GCy +EC,, (21) 


14. All displacements considered are small in magnitude, 
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Eq. (21) includes one more term containing the product of inertia Ing as com- 
pared with Eq. (1). Furthermore, Eqs. (22) and (4'') are essentially ‘the same 


2 
except for the constant coefficient of the ($’)dx term. 


Attaching the column to the skin provides certain restraint, namely 


owt 0 


+ 


Eq. (23) may be used to eliminate one unknown displacement, w. Doing so, we 
obtain 


2 


+(ECyt 4y ET) 


- Ing? 


13 18 


Vv = Vo A Vv 


Xe 


Substituting 


into (24), (25) and denoting 


: 
a 
or 
4 
7 
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Adoption of the same reasoning as for the case of free buckling and application 
of (10) lead the simultaneous equations: 


ty 
+(EC,, Ax - x 4x] 
The amplitudes of deflection, vg, oo" become indefinite when 
fac, dx + 
(28) 
(x,)"da | 


(a) 


dx 
| 
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Considering (E Cw)m = ECw+y 2 © I, = modified warping rigidity, and us- 
ing the notations (12), we may transform (28) into 


(29) 
=Q --- 
Upon expansion, we get 
(299 
where 


In most instances X, = X¢g and Myq¢ = 1, then 


q 
A 
+h, A ff 4), 


If the axis 7 or ¢ is an axis of symmetry, then Ine = 0, and Eq. (29') reduces 
to 


(30) 


(31) 


Eq. (31) is identical to Eq. (19), except J,,,' is used instead of J,, to compute 
the coupling factor? . 


a 
Es 
Ve ~ 
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In application of (29') to (31) care must also be exercised that Py is calcu- 
lated from the modified (E Cy)m and Jm'. 


Column with Known Axis of Rotation 


An example of this category is shown in Fig. 5. As a corner column at- 
. tached to skin in two perpendicular directions, the column is enforced to ro- 


tate about the axis M during buckling. This restriction provides the following 
relations: 


6) 


6) 


(32) 


Applying the similar procedure to that illustrated in the previous paragraph, 
we shall conclude that the uncoupled torsional buckling formula 

G Cul +(E Xe) dr 


P= 
(33) 


may still be used to determine the critical load, provided the following modi- 
fied constants are employed: 


(EC,).= E Ty zm + Iy 

Ky ‘kg 


As seen in Fig. 5, 7 and ¢ are the centroidal principal axes throughout the 
discussion of this paragraph. 


Experimental Evidences 


A series of experiments has been conducted by A. S. Niles for columns 
: with simple supports at both ends. The results were published in Ref. (2), and 
Ps ‘ they agreed closely to the theoretical values. It is reasonable to predict the 
similar agreements for other end conditions. However, further experiments 
pertaining to different boundary conditions with and without skin attachments 
are highly desirable. 
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Formulas for Cy, ‘ty and *, values for some conmonly 


used sections: 


(1) Eqnal-lerred Angle 


(2) I Section 
b 
t,, F 


_ 2ht, + bt, 


4 

b 

| w 12 ht +2bt, 
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(4) Channel Section 


q 
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(3) Z Section 
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(4) C Seetion of Uniform Thickness 


- EY] 
(4° +2¢ *+42be +b?) 
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(5) Hat Seetion of Uniform Thickness 
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DYNAMIC EFFECT OF A MOVING LOAD ON A RIGID FRAME 


R. C. DeHart,! M. ASCE 
(Proc. Paper 1794) 


SYNOPSIS 


A procedure for determining the response of a rigid frame structure sub- 
jected to a distributed load moving at a constant velocity is described in this 
paper. It is demonstrated that a single impact factor is not applicable to all 
parts of the frame. In addition, tables are presented which facilitate the 

handling of a distributed load of variable intensity. 


INTRODUCTION 


In an earlier paper(1) written by the author, the basic dynamic response 
equations for a pin ended rigid frame were developed and an analysis made 
of the response of a rigid frame subjected to time dependent horizontally ap- 
plied loads. This same technique is now utilized to study the case in which 
the load is moving at a constant rate across the top of the frame. Finally a 
comparison is made between the moment caused by a statically applied load 
and a moving load having the same intensity. Tables for simplifying the 
computations required if the integration of equation (8) is to be done 
numerically are also given. 


List of Symbols 


yu" is mass of the column per unit length. 
uY is mass of the beam per unit length. 
L" is column length. 


LY is beam length. 


ete Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
a written request must be filed with the Executive Secretary, ASCE. Paper 1794 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 4, October, 1958. 


1. Structural Research Analyst, Armed Forces Special Weapons Project. 
Washington 25, D. C. 
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Eis modulus of elasticity. 


I is moment of inertia. 
Wj is the vibratory frequency of the frame in the ith mode. 


- sin + D, sinh 


Equations of Motion 


For the loading and frame arrangement displayed in Figure 1 the expres- 
sion for the generalized coordinate is(1) 


v t 
A, sin + B, cos wt + (r) sin w, (t -t) dt 
aw 
i 
(ay)? 
i | i i 


* T, is the kinetic energy in the i” mode, 


BI” 

ET 

For unsymmetrical modes Ee 
u 

i 

v 

u vL u 

v 

vL u u 

For symmetrical modes Sal 
v 
u vL 
i i 
at 
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Figure l. Loading and Frame Arrangement 


. is the quantity obtained when the integral of equation (36), Reference 1, 


is evaluated without the . factor included and is multiplied by 


ky 


, 4 is obtained from equation (35), Reference 1, without the - multiplier 
k 
included. 


i 


v 2 


1 


The generalized force for the unsymmetrical modes is 


st 


v | Vv 
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v 
vi (cos k’x cosh sin k’x 


4 i 


v 
sinh k,x dx 


And for the symmetrical modes is 


st 
Q(t). P | cot we (cos = cosh + sin k*x 
4 


cot k 
+ iz sinh k’x | dx 
vL i 
coth 


The integration of equation (2) for the unsymmetrical modes results in the 
following: 


v 
L 
v tan kY 
- tan (sin - sinh k,8 tT) cosh k,s T= cos k,s 
i 
v 
tan k 
+ +1 (4) 
tanh ks 


Substituting equation (4) in (1) yields 


tan 


v 
cosh k,st 


vL 


vL" 
i 


sin wy (t -t) at 


Upon integration of equation (4), the following is obtained: 


(2) 

Vv 
vL 

tan 

vL 

(3) 

v 

vL 
tan k is 
(5) 
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¥ ar k,8 sin + w, sin k,st 
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wi’ sinh k,st k,8 sin w,t an 
iz 
4 
v (6) 
tan k (= k, st cos 
vL_ Vv 
(xo) +4 7 
cos k, st cos w,t tan l<cos wt d 
— Similarly for symmetrical modes: 
v v 
v w, Sin k,st - k,s sin 
vi" w, Sinh k,st - k,s sin 
coth (k¥s) + 
cos = cos u,t cot lecos w,t 
vi 
pa Equations (6) and (7) are applicable only for the time required for the load Be 
:. to traverse the beam. This time period will be referred to as the first era. a 
5 . . During the second era, the entire beam is loaded and the expression for ay for - 
symmetrical modes becomes 
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qy A, sin + B, cos (8) 


For the unsymmetrical modes equation (8) applies except Q; is zero dur- 
ing the second era. 


Values for q and Qa’ are obtained at the end of the first era. This infor- 


i 
mation permits the determination of the constants A; and B;. 


TABLE 2 
Values of A;, By, and Time to 
Second Era Meximums and Minimums 


Second " Not determined since moment contribution is es 


deci of 
ntic vr 
entialily nile 


| 

| 

First Unsyrmetrical] + 0.2172 - 0.5698 0.576 seconds | 
| 


Third " + 9.001411 + 0.0097239 0.34 

Fourth 0.0005064 + 0.000L9LE 

#* For unsymmetrical modes, time to second era maximum inimum qi i ™ ) 
is time to maximm second era negative and positive moment, respect ye 


These constants for the various modes are displayed in Table 2 along with 
the time at which the initial second era maximum and minimum qj values oc - 


cur. Later maximum and minimum q; values would be the same since damp- 
ing has not been considered. 


Determination of Bending Moment 


The following relation between bending moment Mj, generalized coordinate 
qj and mode shape yy will be utilized. 


= 
d 
; 
Time to Initial ' Time to Initial = 
Mode Ay B. Second Era, Max- Second Era, Min- d 
imum qi (& Moment)} imum q; (& Moment) 
First Symmetrical 0.2805 = 022615 | 2332 Seconds 0.432 Seconds 
Second " + 9.00125 | = 9.903024 0.322 | 0.340 
‘ 
4 
a 
a 


I is the moment of 


Figure 2e 


v 
= 1,93 


vLY 
3-16 


v 
VL 
3-98 


Vv 
7604 


v 
vL 
10,1 


where E is the modulus of elasticity 
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Ox 


inertia 


Bent Dimensions and Load Intensity 


Vv 


0,333 


| 360" | 


lb. sec” 


in 


= 39h in? 


For the unsymmetrical modes 


1” = 1043 in 


1.05 


= 3.22 
4.07 


u 
= 7.18 


u 
kel 10.3 


Velocity = 1200 in/sec 


The coordinate system and loading are displayed in Figure 2. The bent 
under consideration is the one treated in Reference 1 from which most of the 
basic information listed below is extracted. 


w. * 9,98 rad/sec 


= 94.1 rad/sec 


= 150 rad/sec 


L468 rad/sec 


969 rad/sec 


| 


| 
20 1b/in 
Load intensity 
q 
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For the symmetrical modes 


1.87 T = 1.12 31.8 rad/sec 

= 1.83 = 1.87 1 

= 3.59 kL" = 3.66 "77.0 w* 122 rad/sec 
= 5.07 T, = 3013 = 234 rad/sec 


= 6,82 422 rad/sec 


a 
A 
5 
z 


Symmetrical | 


0.3 
Time in seconds 


Figure 3. Time Dependent Moment for Symmetrical Modes--Left End of Beam 


Time dependent moment for the left end of the beam for the first three 
symmetrical modes and the first, third, and fourth unsymmetrical modes are 
given in Figures 3 and 4, respectively. The fourth symmetrical, and second 
and fifth unsymmetrical modes were, for all practical purposes, not excited. 
It is apparent that bending moment contributed by higher modes oscillates 
rapidly relative to the oscillation of bending moment arising from the lower 
modes. A change in a detail such as the amount of end restraint at the column 
footing produces a greater change in the frequency of one mode than another. 
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Figur me Dependent Moment for Unsymmetrical Modes--Left End of Beam 
Figure 4. Time Dependent } y 


For example, the influence of changing from the pin ended frame used as an 
example here, to a fixed ended frame of identical proportions, is given in 
Table 1. On the other hand, small changes in the frequency of a mode make 
only small changes in the moment contributed by the mode. 

It, therefore, appears to be impossible, in the practical case, to determine 
the exact phasing of the modes. This immediately leads to the concept that, 
except for perhaps the first symmetrical and unsymmetrical modes, the 
maximum positive and negative moment from each mode can occur simul- 
taneously. This procedure would give the upper bound for both positive and 
negative values of moment. For the five unsymmetrical and four symmetri- 
cal modes examined, the upper bound for moment at the left end, the center, 
and the right end of the beam would be 1,700,000 in. lb, 2,370,000 in lb, and 
1,700,000 in. lb, respectively. All of these upper bounds occur in the second 
time era. 

Comparable maximum values for the case in which a static load of 240 
pounds per linear inch is applied on the full length of the beam, are 1,500,000 
in. lb, 2,370,000 in. lb, and 1,500,000 in. lb. If more modes were used, it is 
obvious that the center moment, under the moving load conditions, would 
probably be somewhat more than the static center momeni. 

The conclusion may now be drawn for the bent studied herein that, at the 
center of the beam, the dynamic bending moment arising from a uniformly 
distributed load moving at constant velocity is, for all practical purposes, 
identical to the maximum static bending moment produced at the center of 


& af 
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TABLE 1 


Influence of End Conditions on 
Frequency of One Bay Single Story Bent 


Mode Pin End Freouency Fixed End Freauenc} 


31.8 RAD/SEC 33.7 RAD/SEC 


First Symmetrical 


Second " 170. 
Third " 25. 5 
Fourth " 08. 20 


First Unsymmetrical 


Second " 


the beam by a uniform load having the same intensity. This is the same con- 
clusion that Inglis(2) came to in his study of a simple beam subjected to the 
loading conditions described above. 

On the other hand, under these dynamic loading conditions, the increase 
over static bending moment at the ends of the beam is approximately 11%. 
This obviously indicates that impact factors, if they are to be used to account 
for the difference between dynamic loading and static loading, will vary for 
different portions of the structure under examination here. 

Of further interest is the fact that, had sidesway been eliminated, excita- 
tion of the first unsymmetrical mode would be greatly reduced, and the left 
and right end upper bound bending moments would be nearly the same as for 
the static load case. One would, therefore, conclude that in this instance the 
moving load produces a bending moment essentially no different than the 
maximum static bending moment arising from the application of a uniformly 
distributed load to the beam. 

It should not be overlooked, however, that the velocity with which the load 
traverses the beam relative to the frequency of the various modes of vibra- 
tion, particularly the first symmetrical and unsymmetrical modes, deter- 
mines the magnitude of the dynamic response. 


Numerical Integration of Equations of Motion 


Attention is now given to the development of information necessary if the 
equations of motion are to be integrated by numerical means. It is convenient 


to use the Lagrange equation of motion written for the ith mode, in the follow- 
ing form: 


6 
Third " 150. 200. 33 — 
Fourth 373. 391. 5 
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(9) 


2 % a G; is usually referred to as the generalized 
1 


The quantity 


mass for the it mode and will be denoted by M,. The quantity 2U (kY) er” 
may be thought of as corresponding to the spring constant for a single degree 
of freedom system, and will be denoted by Ky. Obviously, no particular 
terminology is necessary for these factors. 

Since U; and T; have equal values, all quantities, with the exception of Q;, 
are obtainable. 

Q; is defined by equations (2) and (3). For an arbitrary distributed load 
(time dependent and space dependent) Qy can be plotted as a function of time 
utilizing the values given in Table 3. Tsien values are for / yY dz for vari- 
ous lengths of the beam where y’Y is the mode shape equation for the beam. A 


typical Qy versus time curve appears in Figure 5. Equation (8) is now writ- 
ten, using the terminology given above, 


(10) 


Equation (10) can be integrated by numerical or graphical means and 
values of the moment contributed by the ith mode determined from the rela- 
tion between bending moment and mode shape. 

If the load is of constant intensity, there appears to be no advantage in 
numerical integration. However, if the load intensity varies, for example as 
given by Figure 6, numerical integration would be an attractive means of ac- 
complishing the integration of the equations of motion. 


Load Intensity 
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SYNOPSIS 


Closed form solutions are obtained for the deflections, residual deflections, 
residual membrane tensions, and other quantities of engineering interest for 
infinitely long clamped rectangular plates with large deflections, under uni- 
form pressure. The analysis assumes infinite rigidity, in the plane of the 


plate, of the boundary supports. Charts are provided for convenience in 
solving the equations. 


INTRODUCTION 


The so-called plastic design of beams and frames has received consider- 
able attention in recent years. In extending these methods to plate analysis 
several difficulties arise, not the least of which is that the small deflection 
theory of plates due to Lagrange(1) yields a limit strength much lower than 
that actually developed by plates of practical proportions. This is because, 
as the plate continues to deflect, membrane tensions arise in the middle plane 
of the plate which add to its load-carrying capacity. A good estimate of the 
ultimate loads on plates, neglecting membrane action, can be obtained by the 
rupture-line theory of Johansen.(2,3) While it may be permissible to neglect 
membrane action in certain practical applications, all plates develop mem- 
brane forces for large enough deflections. Any realistic plastic analysis 
must, therefore, take into account the effect of membrane tensions in the 
middle plane of the plate. 

The basic equations for plates in the large deflection region have been de- 
rived by von Karman.(1) An exact solution of these quations in the elastic 
case was obtained by Levy.(4) When one attempts to make an elasto-plastic 


Note: Discussion open until March i, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1822 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 4, October, 1958. 

1. Southwest Research Inst., San Antonio, Tex. 
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analysis of rectangular plates using the von Karman equations as the starting 
point, the mathematical difficulties are prohibitive. One serious difficulty is 
that rectangular plates with clamped edges yield at the middle points of the 
edges, and then the plasticity spreads gradually to the corners and to the 
center of the plate. Thus, at any given stage in the analysis, the boundaries 
may be partly elastic and partly plastic. Thus, in addition to the inherent 
mathematical difficulties of von Karman’s non-linear equations, one has to 
consider the complex situation of mixed boundary conditions. 

The particular difficulty associated with mixed boundary conditions can be 
avoided if the plate is assumed to be infinitely long. It is then possible to 
consider an elementary strip of unit width (Fig. 1) and length) , and 
von Karman’s equations reduce to a single linear differential equation:(1) 


ax4 D v d x¢ 


Eh3 
where D = 12(1-v2) is flexural stiffness of the plate, q is the load on the 


plate (assumed uniform in the present analysis) and N, is the membrane 
tension which, in this case, is constant throughout the plate. As one is con- 
sidering an elementary strip, the boundaries may be assumed to yield through- 
out at the same time and the complexity of mixed boundary conditions thereby 
avoided. Furthermore, the plastic-hinge concept, so effectively utilized in the 
case of beams, may be readily extended to the plate problem. 

Clarkson(5,6) has solved the above equation for the elasto-plastic range 
by a power series method, in which he uses two terms of the series to obtain 
the final solution. In the present paper, the basic concepts of Clarkson are 
adopted, but a closed form and therefore “exact” solution is obtained, within 
the limits of the assumptions stated below: 


1. The stress-strain law. The relation between stress and strain for first 
loading, unloading, and reloading is given in Figure 2. The analysis may 
therefore be expected to apply to materials which exhibit little or no strain 
hardening, or have a large range of pure yielding. 
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2. The yield condition. The Hencky-Mises yield condition(7) is adopted in the 
analysis. For a biaxial stress condition this becomes 


2 2 
Gx2 + - Ox Fy + 30 xy o>; (2) 


where @, is the yield stress in pure tension. For the infinitely long clamped 
plate consideredhere ¢ 
y x xy 


(see page 2 of Reference 2). Equation (2) thus simplifies to 


os 
Sx (3) 
3. Moment-curvature relation. The moment-curvature relationship is simi- 
lar in form to the stress-strain law. This involves the assumption that a 
cross section of the plate behaves elastically until it has yielded right 
through, that is, until the distribution of stresses is as shown in Fig. 3. 
From there on its resisting moment remains constant at and the sec- 
tion forms a plastic hinge exactly as assumed in the plastic theory of 
beams. The plate is free to rotate about the hinge with no resulting change 
in stress at that section. Mp will be called the “fully plastic moment” and 
differs from its value in the beam theory only insofar as it is modified by 
the presence of a membrane tension, and Poisson’s ratio. 


Referring to Fig. 3, one has 


Gs 


Ji-v+v" 


eliminating h, between these two equations 


= 6s ne lew+v* (4) 
2.2 
4 [i-v 6; “h 


The first plastic hinges form at the edges of the plate, where the bending 
moment is the maximum. As the load is increased, a plastic hinge forms 
along the center line of the plate. The load corresponding to the formation 
of a center-line plastic hinge is called the “limit load.” (Actually the plate 
will continue to carry load beyond the “limit load,” until the whole width of 
the plate is plastified. But this phase is beyond the scope of the present 
paper). 


4. Other assumptions are the same as in large-deflection plate theory in the 
elastic range. Briefly, the material is supposed to be homogeneous, 
isotropic and to obey Hooke’s law until it yields. The deflections may be 
many times the thickness of the plate but small compared to its lateral 
dimensions. 
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On integrating equation (1) twice with respect to x, one has 
2 
dx 2 2D D 2 


where cy and Cg are constants of integration. Referring to Fig. 1, and mak- 
ing use of the relationship, 


the constants of integration are easily determined and the above equation be- 
comes 


é 


— + qx + Mo 5 
axe D 2 2 D (5) 


The general solution of this equation is (see page 12, Reference 2) 


22 
sinh , cosh 2ux_, gi>x 
8u &u “D 


(6) 
la 


from which 


3 2 
2u_ cosh _2ux 4, Co _2u_ sinh 2ux ge =((7) 


where 
42 
Nxt 

(8) 


u being the “tension parameter” and c, and cg are arbitrary constants. Sil 
The boundary conditions, valid in both the elastic and plastic ranges are: 


w=0 at x=0O andx=t 


= =o at x= Y2 


On substituting these conditions and then transferring the origin of co- 
ordinates to the center of the plate (Fig. 4) one finds: 


cosh 2ux 22 
cosh u 32u°D Buep 


dw . sinh 2ux gtx (10) 


ax cosh u 
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h= thickness of plate 


A q B 
w 


Compression | Tension 


Fig. 3 Fig. 
where 2 
co = + M 
1 
2 


The condition for the calculation of »< = 


is obtained by considering 


the change in the distance between the supports A and B (Fig. 4). Assuming 


that this change is zero, i.e., the edges are rigid enough to present any “pull 
in,” this condition may be started as follows: 


1 L 
2 2 ix 2 
© 


which becomes 
Eh ox 
E! ax 
On substituting the value of N, from equation (8), this reduces to 


(— >) ax (11) 


On page 12 of Reference 2, it is shown that for the elastic case where the 
boundary slope is zero, M, has the value 


coth u (12) 
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When a plastic hinge is formed at the edges, the end moment is the fully plas- 
tic moment given by equation (4). Mp may be expressed in terms of the 
“plate strength parameter” B introduced by Clarkson and defined as follows: 


Bs 12 (13) 
1» + E ne 
Mp may then be written 
| = _3hD Be 14 
P 12 9B (14) iis 
Thus, after the formation of plastic hinges at the edges, the constant Co -_ 


becomes 


(15) 


Substituting (15) and (10) into (11) and integrating, the condition for no “pull 
in” may be written in the following form: 


9B 1 d; 9B2 ub 6 


the “non-dimensional load” (17) 


Equation (16) is a quadratic for P for given values of u and B. Figures 
5 and 6 give a plot of P against B for various values of u. In Fig. 5, P is 
plotted on a logarithmic scale. Table I contains the same information as 
given in Figures 5 and 6. 

The deflections after the formation of edge plastic hinges may be obtained 


by substituting the value of cg from equation (15) into equation (9). Thus one 
finds: 


ccsh 2ux ot 


The maximum deflection may be obtained by letting x = 0 in (18). It is con- 
venient to express this in non-dimensional form as follows: 


(i-cosh u 2 »+3/u (i uly (19) 
u-coshu 9B u~cosh u 


3h 
u 
16u'D 
et 
(16) 
in which = 
@ « . — 1... 4, 
u? 6u ° lu°cosh*u 
tanh u 1 1 
2u 3u ou “cosh u 
3 2 
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If corresponding values of B, u, and P are substituted in (19), the maximum 


w 
deflection may be obtained. In Table I are shown values of a this calcu- 
lated. 


An inspection of Table I shows that, for all practical purposes, the non- 


w 
dimensional deflection aa is independent of P and B, and may be con- 


w 
sidered as a function of u alone. When average values of > are plotted 


against u (Fig. 7), it is seen that the linear relation 


(20) 


w 
closely expresses the relation between u and 
The end slopes corresponding to any load after plastic hinges have formed 
at the edges can be obtained similarly from equation (10). It is convenient to 


express the value of the end slope in “non-dimensional” form thus: 


1x n By‘ 


Fig. 8 gives a plot of @ against P for various values of u. P is plotted ona 
logarithmic scale. 


4+ ont 


(21) 


Onset of Edge Plasticity 


The load when plastic hinges form at the edges may be determined from 
the condition that the elastic fixed-end moment equals the fully plastic 
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TABLS I (LOADING) 


u B P Wn/h A emarks 
1 1.5), 82.6 0.367 Center Plestic Hinge 
156.2 2369 Edge Plastic Hinge 
2368 a 
2 3.25 22361 e747 1.713 Center Plastic Hinge 
5.00 282.1 0756 1.208 
9.14 396.9 0739 0 Edge Plastic Hince 
a 
3 5.h6 505.8 1.108 20765 Center Plastic Hinge 
10.00 648.5 1.118 1.491 
15.00 784.5 1.098 0.118 
15.42 804.2 1.107 0 Edge Plastic Hinge 
1.108 a 
4 8.21 981. 1.472 3.972 Center Plastic Hinge 
10.00 1050.5 1.1482 3-475 
15.00 1213.9 1.494 22203 
20000 1354.6 1.488 0.977 
23095 1455.4 1.472 ) Edge Plastic Hinge 
1.482 a 
5 11.60 1730.5 1.826 5.238 Center Plastic Hinge 
15.00 1867.8 1.6h4 4.361 
20.00 2037.0 1.856 32190 
25200 2185.0 1.858 2.074 
30.00 2318.9 1.851 0.971 
34.38 2426.6 1.836 9) Edge Plestic Hinge 
1.847 a 
6 15.66 2819.8 20179 6.543 Center Plastic Hinge 
20.00 3004.6 20199 
25200 3184.0 2.212 4.403 
26230 3226.9 2.215 4.131 
30.00 3342.4 2.219 30371 
40.00 3618.9 20213 1.372 
46.89 3786.1 20197 Edge Plastic Hinge 
20205 a 
7 20.40 4317.6 2.530 7.873 Center Plastic Hinge 
25000 4525.1 2-59 6.824 
26.30 4577.8 2.554 6.548 
30.00 4717.8 256k 5-793 
40.00 5045.0 20577 3.889 
50.00 532h.2 2.575 2.071 
61.52 5603.6 20555 0 Edge Plastic Hinge 
2.554 a 


a@- average for constant u 
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60-00 


90000 219 10.704 

90.00 | 9.37 
100.00 8.0L7 
120.00 606.8 
140.00 3° 32127 
160.00 06745 


a = averare for constant u 
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3 
4 sADLE (LOADING) - ) 
4 u B P Wn/h A Remarks J 
8 9.217 Center Plastic Hinge 
— 9 «106 
8.313 
6-402 
4.660 
i 22990 
0 Edge Plastic Hinge 
a 
9 31.91 £€10.0 10.575 
‘ 40.00 8.981 
50.00 
©0200 54626 
< 200 2 62 
90.00 1.02 x 
a 
10 30.0 35.9 3-5€2 11.939 = 
40-00 109 3.566 11.678 
£00 179 60 3613 9.872 
60.0% 19067 3-631 
79.00 3-641 2677 
90.00 OS el 3-646 3.883 
100.00 1065 26 3-641 2-491 
110.90 130500 30632 1.109 
™ 115.05 440563 3.621 0 Edge Plastic Hinge 
3.62, a 
12 2030404 14.688 Center Plastic Hinge 
2061326 13.703 
A 
al 3 
Edge Plastic Hinge 
4.331 a 
q 
ae 
a 


(20) (28) 
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moment. Thus, from equations (12 and 14): 


1822-13 


12 ue tanh u ) 
te 
u= tanh u 9B 


where 


a Dh 


The relation between B and u when edge plastic hinges form can be derived 
by substituting the value of P from (22) into (16). In this way one finds: 


27 anh@u tanh u 3 


u-ta 


8 2 
utah» | @, utanh’u é 2 ri} 


u-tanh u) (u-tanh u)2 8 


ylé tanheu D2 u-tanh u ui ® _ 0 
12 (u-tanh u)2 4 v=tanh u 48 (23) 


where $4, $2, and @g, are the same as in equation (16). 

Equation (23) gives related values of B and u which must be substituted 
in equation (22) to obtain the non-dimensional load P corresponding to edge 
plasticity. Figs. 5 and 6 are shown @ plot of equation (22). Table I gives the 
values from which the curves were drawn. 

It may be noted here that when there is no tension in the middle plane of 
the plate, i.e., u = 0, equation (23) no longer applies, and B is independent 
of u. On setting u = 0 in equation (22), there follows 


P, = 36B (24) 


This is the non-dimensional load for the formation of edge plastic hinges 
if the tension in the middle plane is neglected as in classical plate (or beam) 
theory. Fig. 6 shows a graph of equation (24) and clearly indicates the great 
increase in the stiffness of the plate due to membrane tension. 


Onset of Center Plasticity 


The load of P corresponding to the formation of a center-line plastic 
hinge may be obtained by equating the bending moment at the center to the 
&s fully plastic moment. Differentiating equation (18) and letting x = 0, there 


follows: 
2 
Ix = 6 luScosh u Im? tcosh u 


Setting this equalto Mp = (B - 


) and simplifying, there 


finally results: 


3 
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quc cothu 3Dh ut 
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ud cosh utl 
(25) 


cosh u- 


P, = 12u* (B - 


This is the “limit load” on the plate as here defined. 
The relation between B and u when a center-line plastic hinge forms may 


be obtained by substituting P from (25) into (16). The result may be written 
as a quadratic equation for B2: 


(cosh u-1) u2 (cosh u-1) 
® ue (cosh utl) 3 sh #1)2 
2 + = {cosh + 3 gq, +1 
(cosh u-1) 2 (cosh u-1) 8 
1 


+ , u (cosh ust) , 


12 (cosh u-1)¢ 


(cosh u-1) 8 (26) 


Figs. 5 and 6 show a plot of equation (25), corresponding values of B and 
u being obtained from equation (26). 

It may be noted that when there is no tension in the middle plane of the 
plate, equation (26) no longer applies. Setting u = 0 in equation (25) yields: 


P,, = 48B (27) 


This is the non-dimensional load for the formation of a center-line plastic 
hinge when the tension in the middle plane is neglected. Fig. 6 shows a graph 
of equation (27). 

The non-dimensional edge slope when a plastic hinge forms at the center- 
line may be obtained by substituting the value of P from equation (25) into 
equation (21). One finds then: 


A 


sh u-2 sinh 
( cosh u-2 sin (28) 


2 ut 
Fig. (8) shows a plot of equation (28) 


Unloading 


When unloading takes place after the formation of plastic hinges at the 
edges, the end rotations remain constant at the maximum value reached dur- 
ing first loading. Otherwise the beam behaves elastically. Thus the maximum 
end rotation reached prior to unloading provides the additional boundary con- 
dition required for a unique solution. 


It is evident that the basic differential equation for unloading is (5). The 
boundary conditions (Fig. 4): 


at x=0 


dw 1 
a =o atx 


are still valid. If these are substituted in the solution (6) and the origin 
transferred to the center line (Fig. 4) one obtains equations (9) and (10). 
The value of the constant Co is the same as before, and Mo is now deter- 
mined from the condition 


* 
BY 
2 ue 
B 
-B 
= 
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° tanh u 7 2uDh 
The notation P is used to distinguish the load reached during the unloading 
. process from the load P reached during the loading process. The deflection 


during any stage of the unloading may be obtained by substituting (29) into (9). 
One thus finds: 


cosh 2ux 
Ph ah Ph Px<h 


It is to be noted that a has to be substituted from equation (21) and that the 
value of u to be used in determining @ corresponds to P and not P. The 
maximum non-dimensional deflection can be obtained from (30) by setting 


x =0. 
Wr 7 sh a } c 
I l-cosh u +1 l-cosh u 
h l6u* { u sinh u | z u sinh u ( ) 


The value of u in equations (29) and (31) is obtained by using the condition of 
no “pull in” as given in equations (11). One thus finds: 


Fig. 9 gives a plot of P against o for various values of u. Table II gives 


w 
the same information. This table also gives the values of — calculated 


w 
from equation (31). An inspection of these figures shows that a may, for 


all practical purposes, be _— as a function of u alone. Fig. 10 shows 


a plot of average values “a > for a constant u, against u. It is seen that 


the relation is almost exactly linear and may be expressed by the equation 


Ww 
— = 0.361 u (33) 


It is a remarkable and unexpected result that equations (20) and (33) are 
identical. 


a 
4, 
Substituting from equation (10), one finds 
| 
a 
4 
a 
52 9 3 3 1 } 
= 2 
Bu léu-tanh u 1éu4sinhéu loa 3tanh u lu*sinh@u 
4 
q 
(32) 
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TABLE II (UNLOADING) 


OF 
. 
38 
Pa 
*e 
I 
@ 


3.972 981.9 1.472 Center Plastic Hinge 
3200 1.469 

2.00 1238.9 1.49) 

1.00 352.8 1.489 

0.00 1456.2 1.73 Edge Plastic Hinge 


J 


173007 

1768.7 le 
1921.4 le 
206361 le 
2194.4 
2315.5 1.8 
2426.6 


LW lw lo lw PO 


ASCE 
1 82.6 Center Plastic Hinge 
0.00 156.2 2370 Edge Plastic Hinge 
0369 a |. 
2 1.743 22361 Center Plastic Hinge 
r 
1.00 302.6 0754 
Fe 0.00 396.9 0139 Edge Plastic Hinge 
e750 a 
3 26765 505.8 Center Plastic Hinge 
2.00 591.6 1.12 
1.12 a 
4 
1.48, a 
5 5.238 6 Center Plastic Hinge 
5200 2 
3200 7 
Q 
2.00 
1.00 1 
0.00 6 Edge Plastic Hinge 
1.844 a 
6 6.543 819.7 2.178 Center Plastic Hinge 
6.0 917.0 22190 
520 087 20206 
ia 2u7 e2 2-216 
300 39667 2.219 
4 2.0 536.3 2.217 
0.0 786.1 2.196 Edge Plastic Hinge 
2-204 a 
a = average for constant u 
q 
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TABLE II (UNLOADING) - Cont'd 


Wm/h Remarks 


lee) 


2.529 Center Plestic Hinre 
225146 
22561 
2571 
20577 
2.578 
20574 
22567 
5603.6 22555 Edge Plastic Hinge 
2.562 


ann 
. 


888888ss 


~ 


6291.5 2.880 Center Plastic Hinge 

6558.8 22902 

6766.7 22916 

696.4 26926 

715202 20932 

733003 20935 

7498.7 20934 

7657 07 2.930 

7807 «2 2.922 

7947 ol 2-911 Edge Plastic Hinge 
2.916 


8809.0 30231 Center Plastic Hinge 
8950.4 3.241 
9188.0 30257 
9415.4 3-270 
963267 3279 
980.1 3.286 
10029 «7 32288 
10226.1 30291 
1040.9 3.289 
1057763 32285 
107346 3-277 
10885.5 32267 Edge Plastic Hinge 
3-272 


10 
10 


888834 


11938.9 3.58 Center Plastic Hinge 
12189.7 3-597 

12))6.8 3.612 

1269347 32623 

129307 3.633 

13158.0 3640 

1337507 3.644 

1358329 32616 

13782.8 3.646 

1397223 3.643 

1152.6 3.638 

1323.8 3.631 

3-622 Edge Plastic Hinge 

3.627 a 


88 


Wan 
eeee*# 


8888888888 


a = average for constant u 
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u x“ P 
he 
36 
2 
he 
9 
} 
6 
5 F 
3 
2 
9 
9 
8 
7 
6 
5 
3 
1 
0 
10 


a 


) 
203044 


14.00 20521.0 
13.00 208275 4.308 4 
12.00 2112.0 4.320 
11.00 21110.6 e330 
10.00 21687 .6 4.339 
9200 2195409 
8.00 2221269 4e350 
7200 22461.) 
2270026 4.355 j 
229 30.6 


} 


2394309 


a = average for constant u 


Unloaded Plate 


Certain important quantities when the plate is completely unloaded may be 
obtained from the above equations by letting P = 0. oa 

The maximum residual deflection is obtained by letting P = 0 in (31). 
There follows: 


(in XA(cosh u-l) 


>? residual 2u sinh u (34) 


In equation (34) a is the maximum non-dimensional edge slope reached prior 
to unloading and u is the residual value of the tension parameter after com- 
plete unloading. The relation between these two quantities is obtained by 


setting P = 0 in equation (32). Thus: 


tanh u sinh*u (35) 
3 sinh@u-u tanh u 


Equations (34) and (35) are plotted in Fig. (11). Table I shows the 
figures from which these curves have been drawn. 


Example 


In Figure 12 are shown two curves taken from Mr. Clarkson’s paper.(8) 
One is a plot of non-dimensional deflection against non-dimensional load and 
the other non-dimensional residual deflection against non-dimensional load. 
The curves are for a plate of strength parameter 26.3. On the same chart 
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TABLE III (UNLOADED PLATE) 


RESIDUAL DEFLECTION & TENSION 


are shown circled points representing results obtained by interpolation from 
curves given in the present paper. (Some of these values were also calcu- 
lated and are shown in Tables I and II). The agreement between the two solu- 
tions is clearly close. 

To obtain, for example, the load corresponding to the formation of a plas- 
tic hinge at the center line, Figure 5 is used. For B = 26.3, it is seen that 
Py, = 6400 and u = 8.1. From Figure 7, the non-dimensional deflection is 
found to be 2.95. For permanent set information read in Figure 5, P = 4600 
for u = 7. Entering Figure 8 ‘with these values, a = 6.6. With o = 6.6, 
Figure 11 gives the residual non-dimensional deflection as 0.95. 


CONCLUSIONS 


The question naturally arises as to the applicability of these results for 
plates with length-width ratios other than infinity. Levy and Greenman(8) 
have made a study of this problem for the elastic case by using the exact 
solution of von Karman’s equations for a clamped rectangular plate with 
length-width ratio of 1.5. Their conclusion is that, for the elastic range, 
the maximum stresses and center deflection for a plate of such proportions 
differ less than 3 percent from the corresponding values for infinitely long 


Ce 
> 
CE 
ig 
ve 
a 
- +8 
A eached ig 
Peiar + Tin} D 
© Unloading Residual u Residual ¥r 
- 
ny 
5 377 
2 
30471 200 661 
lar 
265 797 
Vevi 320 
piste) 
* 12.92 5.0 1.2 
4 
C20 
"7 
a 
\ 
Sak 
‘ 
= 


October, 1958 


+t 


dial 


Trt] 
44 


J 
1822-22 EM 4 
* 
passes 
+144 4444 +444 +444 ener :3 ‘ 
= 
- 
- 
+ +444 +444 4444+ $444-44444 + +444 
tH Tt +444 44444444 
3 


ASCE ELASTO-PLASTIC PLATES 1822-23 


plates. While this does not give quantitative figures for the elasto-plastic 
case, it would seem a safe assumption that for design purposes, a plate of 
length-width ration > 3 may be assumed to be infinitely long. Clarkson(5) 

has come to a similar conclusion. 

The formulas derived in this paper are entirely theoretical and based on 
several idealizations which cannot be exactly duplicated in practice. In 
particular, the boundary conditions (whether, for example, the edges can be 
assumed to be infinitely rigid in the plane of the plate) will have an effect on 
the results. It is to be noted that while infinitely long plates will not develop 
membrane stress if the condition of “no pull in” is replaced by one of free 
edge displacement, a plate of finite aspect ratio always develops large mem- 
brane stress whether or not the supports are immovable. Furthermore, in 
many practical applications, such as in ship bulkheads, where a panel of plat- 
ing may be regarded as being in an infinite field of such panels, the rigidity of 
the edges is effectively infinite from considerations of symmetry. The 
present solution may be expected to yield reasonably good results for such 
cases. Solutions for many other boundary conditions can evidently be ob- 
tained by the method presented in this paper. 
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B 1. p? os "Plate Strength Parameter*® 
1-v+v E he 
D - En? Plate Stiffness 
12(1-»*) 
E - Modulus of Elasticity 
Nx Membrane Tension 
Dh Eh 
PE - Non-dimensional Load at Formation of Plastic Hinges at Edges 
Non-dimensional Limit Load 
h - Thickness of Plate 
1 - Length of Plate 
q ~ Uniformly Distributed Load per Unit Length 


1 
N 2 2 
u - () Tension Parameter 


w ~ Deflection of Plate 

Wm - Yaximum Deflection of Plate 

Coordinate Axes 

4 "Non-dimensional Edge Slope" 

Poisson's Ratio 

6s Yield Stress in Tension 
Ox - Stress in x Direction 

sy - Stress in y Direction 
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DEMONSTRATIONS OF PLASTIC BEHAVIOR OF STEEL FRAMES# 


Closure by H. M. Nelson, D. T. Wright and J. W. Dolphin 


H, M. NELSON,! D.T. WRIGHT,” J. M. ASCE and J. W. DOLPHIN. ?—The 
writers are very glad to have the contribution on the effects of shear force on 4 
plastic moments from Mr. Sobotka. Particularly so because it shows that the ios 
difficult problem of shear effects has been approached, independently, by the aa 
method that Baker, Horne and Heyman’ also favoured and which the writers ; 
think is the best yet available for practical use. 

However, their equation (1) i.e. 


by Mr. Sobotka’s approach becomes simply 


2 


and they do not see the point of rewriting this in the much more complicated 5 
form of the contributer’s equation (9). I 

The interaction diagram obtained from (1) above is, of course, identical F 
with that in the discussion and it is useful to have this diagram formally 
presented. For design work a different diagram is required for each rolled 
section, and in practice the authors favour equation (1) or the empirical rule 
(3) suggested by the Lehigh team. 

The writers find that in their paper they ascribed this “lower-bound” 
approach to Horne... In fact it appears to have been first published by Heyman 
and Dutton in Welding and Metal Fabrication 1954 and subsequently by fe 
Longbottom and Heyman in the Proceedings of the Institution of Civil a 
Engineers, Part II, August 1956. 

It may be of interest to add the following additional reference: 

The full plastic moments of sections subjected to shear force and axial 2 

load. M.R. Horne, British Weld. Jr. 5 (1958) April. z. 


a. Proc. Paper 1390, October, 1957, by H. M. Nelson, D. T. Wright and J. W. 
Dolphin. 


Lecturer in Civ. Eng., The Univ. of Glasgow, Scotland ; 
. Asst. Prof. of Civ. Eng., Queen’s Univ., Canada | 
Prof. and Head of Dept. of Civ. Eng., Royal Military College, Canada aie | 
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ANALYSIS OF CONTINUOUS BEAMS BY FOURIER SERIES* 7 
Closure by Seng-Lip Lee ; | 
q 
SENG-LIP LEE,’ A.M. ASCE.—Mr. Medwadowski states that the advantage [i 
of having the deflection function and its derivatives continuous over intermed- > 


iate supports given by the Fourier series analysis results also by employing 
Heaviside’s step-function. That the latter is not quite true can be shown by 
means of the following illustration. 

Consider the continuous beam shown in Fig. a. In order to determine the 


maximum value of Ry produced by a moving uniformly distributed load ex- ‘ 
tending over a distance of 1/2 1, Eq. (e) yields 
eo (s+d)£L 
R, = (51 45 Sin 37) sin dx (D-/) 
Mz] 
3 4 (4/ se 45 Sin [ (s +z) ns] = O 
Taking only the first two terms of the series yields “a 
/ - 1411 7§ q 


The solution of which yields S = 0.0639. The corresponding value of R, is 


- 45 sn 24” ) (Cos 6.063900 Cos 


wl 


a. Proc. Paper 1399, October, 1957, by Seng-Lip Lee. 
1. Associate Prof. of Civ. Eng., Northwestern Univ., Evanston, Ill. 
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Fig.a Maximum value of R, due to uniformly distributed load 
extending over a distance of $/ . 


If Heaviside’s step-function were used in this case, Eq. (D-1) would be 
replaced by 


R = 3 (x-£)*+ 


-7[-F XH, (x-$t)*+ | de (D-2) 


523 


Although Eq. (D-2) gives R, in a single expression, the integrand is actually 
discontinuous at x = 1/3 1 and x = 2/3 1. In this problem, it is obvious that 
S = 1/2 <2/3 for which Eq. (D-2) takes the form 


(+H rk G+ 
(4] 
(- $x’ + 35 dx be 
2/4 (- 527+ ex) dx| 


For S + 1/2< 2/3, Eq. (D-2) yields a longer expression for Ry with additional 
integrals. In problems where the limits of integration are not known a priori, 
the choice of one or the other has to be made before differentiating to de- 
termine the value of S corresponding to (R;) max. If the wrong choice were 
made, the procedure has to be repeated. This difficulty is avoided by the use 
of (D-1) where the integral is truly continuous between S¢<x=(5++)f That 


/ 
; 
/ 
J R, | 
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where 
{° 
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the numerical work involved in the solution of the problem under consideration 
by Eq. (D-1) is less than that if Eq. (D-2) were employed is not difficult to see. 
Furthermore, in problems where the limits of integration extend over several 
supports, Eq. (D-2) yields additional integrals whereas the integrand of Eq. 
(D-1) remains unchanged. 


The writer agrees with Mr. Medwadowski that Eq. (15) may be derived si 

from Eq. (11) by the direct substitution of the prescribed deflections of the 

intermediate supports 4(@,,4)-=-d¢,. The reason for using Eq. (13) with the fs 


aid of Castigliano’s theorem can be explained by the example shown in Fig. b. 
To determine the reactions of the springs on the beams loaded as shown, 
the center reaction R, is taken as the redundant quantity. For this case, Eq. Be 
(13) takes the form E 


2 
rs{ Kr 

‘ 2 / 32 2 


where K =k 3 is the spring constant and the values of Ry, the spring 


reactions, are shown in Fig. b in terms of Rj. The application of the theorem 


of minimum energy leads to Bei: 
AU _ 
dR, 
or 
) 
~ / 
— (Psin 3 / 
4 x 2k. ( 2 R, P) = oO 


0.688 + 24.0/g 
| + 72.0/g 


/ 


For infinitely stiff springs, k --c° for which Ry = 0.688P whereas for very 
flexible springs, k-+0 for which = 


It should be observed in this case that the energy method makes it un- 
necessary to consider the deflections at the spring supports and lends itself S 
to a simpler solution, whereas the other approach involves a somewhat length- ‘7 
ier procedure, The above example was a part of the original draft which was ee 
subsequently shortened to the form presented. id 


Fig.b Continuous beam on spring supports. 
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Fig.a Maximum value of due to uniformly distributed load 
extending over a distance of #/ . 


If Heaviside’s step-function were used in this case, Eq. (D-1) would be 
replaced by 


2.3 
{a [- 4, (x $)°+ 


Although Eq. (D-2) gives R, in a single expression, the integrand is actually 
discontinuous at x = 1/3 1 and x = 2/3 1. In this problem, it is obvious that 
S = 1/2 <2/3 for which Eq. (D-2) takes the form 


CH ae 


For S + 1/2< 2/3, Eq. (D-2) yields a longer expression for Ry with additional 
integrals. In problems where the limits of integration are not known a priori, 
the choice of one or the other has to be made before differentiating to de- 
termine the value of S corresponding to (Rj) max. If the wrong choice were 
made, the procedure has to be repeated. This difficulty is avoided by the use 
of (D-1) where the integral is truly continuous between S4<x=(5++)/6 That 
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the numerical work involved in the solution of the problem under consideration 
by Eq. (D-1) is less than that if Eq. (D-2) were employed 1s not difficult to see. 
Furthermore, in problems where the limits of integration extend over several 
supports, Eq. (D-2) yields additional integrals whereas the integrand of Eq. 
(D-1) remains unchanged. 

The writer agrees with Mr. Medwadowski that Eq. (15) may be derived 
from Eq. (11) by the direct re ee of the prescribed deflections of the 
intermediate supports 4(@,4)--¢,,. The reason for using Eq. (13) with the 
aid of Castigliano’s theorem can be explained by the example shown in Fig. b. 

To determine the reactions of the springs on the beams loaded as shown, 
the center reaction R, is taken as the redundant quantity. For this case, Eq. 

. (13) takes the form 


/ 2 3 R? 


where K = k 3 is the spring constant and the values of Ry, the spring 


reactions, are shown in Fig. b in terms of Ry. The application of the theorem 
of minimum energy leads to 


=O 
dR, 


solving which yields 


0.688 + 24.0/g 
/ + 


R 


/ 


For infinitely stiff springs, k -~co for which R; = 0.688P whereas for very 
flexible springs, k-+0 for which R, = *. 


It should be observed in this case that the energy method makes it un- 
necessary to consider the deflections at the spring supports and lends itself 
to a simpler solution, whereas the other approach involves a somewhat length- 
ier procedure. The above example was a part of the original draft which was 
subsequently shortened to the form presented. 


Fig.b Continuous beam on spring supports. 
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Mr. Medwadowski questions the correctness of the statement that a given 
load and the reactions of the intermediate supports can be approximated by 
Fourier series of the type given by Eq. (4). The writer differs with him on 
this point and wishes to point out that the well known technique of representing 
load functions by Fourier series has been used by other authors.2 The 
laborious part in the analysis of continuous beams lies in the determination of 
the deflection and slope functions including the influence lines for reactions. 
To these ends, among others, it has been demonstrated that the Fourier 
series analysis yields simple solutions which converge rapidly. Once the 
redundant reactions are known, the shearing force and bending moment 
functions are best derived by the classical procedure of integrating the load 
functions. 

The writer wishes to thank Mr. Medwadowski for his discussion and 
pertinent questions in regard to the analysis of continuous beams of constant 
flexural rigidity by means of Fourier series. He wishes to express his 
appreciation to Mr. Medwadowski for the opportunity of elaborating on some 
pertinent points of the paper. He is also grateful to him for mentioning the 
paper by Saibel and D'Appolonia. The omission was not intentional. 


2. See, for instance, Bleich, F., “Theorie und Berechung der Eisernen 
Brucken,” Julius Springer, Berlin, 1924, p. 51. 
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SEA BOTTOM PRESSURE FIELDS PRODUCED BY YAWED VESSELS# 
Corrections 


CORRECTIONS: —In the second footnote at the bottom of page 1496-1 and 
continued on bottom of page 1496-2 the phrase that presently reads “at a speed 
(Vgh,” should read “at a speed less than 1/3 Vgh.” The factor b/c appearing 
on the right-hand side of Equation (3b) should be replaced by b/a. The lower 
limit on the integral of equation (3c) should read o instead of A. Also in the 
caption of Figure 2 the ratio a:b:c which now reads “a:b:c = 10.50 = 1:11:1” 
should read “a:b:c = 10.50:1.11:1.” 


a. Proc. Paper 1496, January, 1958, by P. M. Fitzpatrick. 


1. Hydrodynamics Branch, Research Div., U.S. Navy Mine Defense Lab., 
Panama City, Fla. 
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EFFECT OF DEFLECTION ON LATERAL BUCKLING STRENGTH 


1 2 


Corrections by J.W. Clark and A.H. Knoll, 


CORRECTIONS. — 


Page 2, line 2: The reference should be to Pettersson, 
rather than Pattersson. 


Page 11, line following Eq. (17): The variable eliminated 
should be v, notv. 


Page 11: The number for Eq. (18) should follow the 
equation immediately above the one that is now marked 
as Eq. (18). 


Page 12, line 3: Change C to C,. 


. Proc. Paper 1596, April, 1958, by J. W. Clark and A. H. Knoll. 
. Asst. Chf., Eng. Design Div., Alcoa Research Laboratories, New 
Kensington, Pa. 


. Research Engr. Eng. Design Div., Alcoa Research Laboratories, New 
Kensington, Pa. 
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INCREMENTAL COMPRESSION TEST FOR CEMENT RESEARCH 


Discussion by Keith Jones 


KEITH Jones, ! M. ASCE—The civil engineering profession is in urgent 
need of a more complete and detailed understanding of the physical behavior 

of concrete and its constituent materials with particular regard to the elastic 
and creep properties of those materials. That need is accentuated by the in- 
creasing use of thin-shell and prestressed concrete designs. 

The incremental compression test, as devised and described by Mr. 
Hrennikoff, offers a new avenue of approach to the investigation of the 
macroscopic properties of cement and concrete. Further investigations by 
means of this test, particularly if applied to concrete, can be expected to 
yield information applicable to design. It is probable, though, that the 
phenomenon of creep is of a basically molecular nature and therefore fully 
explainable only in terms of molecular movement perhaps related to 
crystalline behavior. 

Previous tests, to the best of this writer’s knowledge, have described 
creep phenomena under long-time sustained loading and the recovery of 
creep deformation has not been so strikingly evident from those tests. Such 
tests do, however, more nearly duplicate the conditions actually experienced 
by concrete in structures. As a more complete understanding of creep and 
related effects such as a variable Poisson’s ratio is gained, criteria may be 
established which would enable the prediction of long-term structural behavior 
of any particular concrete based on findings from short-term incremental 
compression tests similar to those of Mr. Hrennikoff. 

The variation of “ during creep as compared with the constancy of that ratio 
during change of load is, to this writer, the most suprising and perhaps the 
most significant fact presented in the paper. From Fig. 13 of the author’s 
paper it appears that “, measured for a particular period of creep, is related 
to the amount of creep taking place during that period. The top and bottom 
graphs of Fig. 13 show that, in general, the longer abscissas of creep strain 
accompany the higher ordinates of u. In comparison, the constancy of yu 
during change of load is accompanied by substantially equal changes of load 
and of strain. It would be interesting to see a chart similar to Fig. 13 
applied to a test of a dry specimen similar to Fig. 14 or 15 where the creep 
strain for each period is smaller and more nearly constant. 

A second significant fact is revealed by comparing the area between the 
loading and unloading curves of moist specimens with that of dry specimens. 
That hysteresis loop represents a net expenditure of energy which is much 
greater for moist than for dry specimens. It is probable that a complete 


a. Proc. Paper 1604, April 1958, by A. Hrennikoff. 
1. Engr., Commissioner’s Office, Bureau of Reclamation, Denver, Colo. 
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accounting of the energy represented by the hysteresis loop will eventually be 
available to aid the rational and mathematical explanation of creep and its 
attendant phenomena. 

Mr. Hrennikoff is to be congratulated on his initiation of a new test aimed 
at determining the fundamental properties of our most economical, widely 
used, and least understood structural material. His departure from conven- 
tional concepts based on elastic theory where the effects of creep are subor- 
dinate and his substitution of tests in which creep is allowed to assume equal 
importance with elastic deformation may be expected to contribute new and 
more usable and accurate design criteria as those tests are extended in scope, 
usage and acceptance. 
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LIMIT ANALYSIS OF SIMPLY SUPPORTED CIRCULAR SHELL ROOFS” 


Corrections 


CORRECTIONS— Because of errors and omissions which occurred inadvert- 
ently in the printing of Proceedings Paper 1706, a listing of necessary cor- 
rections as given below is considered advisable to permit proper understand- 
ing of the paper: 


On page 1706-4, after the 3rd paragraph ending “and these are pictured in 
Fig. 3”, insert the following: 


Equilibrium.— A precise theory of cylindrical shells must be concerned 
with at least 8 stress resultants, namely 3 membrane forces, 2 radial shears, 
1 twisting moment and 2 bending moments. However as reported by 
Timoshenko (Reference 3, page 447) and in ASCE Manual No. 31 (Reference 
4, page 43), practical results in the elastic range for the type of shell herein 
considered indicate that the only significant stress resultants are the 3 mem- 
brane forces, N,, Ng, S, the transverse radial shear, Q, and the transverse 
bending moment M. We shall here make the assumption that the remaining 
stress resultants may be neglected in the plastic as well as in the elastic 
range. Further, the usual assumption that the radial shear does insignificant 
work is also made herein. Therefore, the radial shear is a reaction, rather 
than a generalized stress (Reference 11); it does not appear in the yield con- 
dition, and it may be eliminated from the working equations. 

With the coordinate system of Fig. 1, the equilibrium equations using 
dimensionless stress resultants are 


(b), 


Cu 


(c). 


Qa 


For the case of the shell with simply supported ends and free edges, the 
stress resultant boundary conditions are noted below. 
(1) Along the free edges g =* @ 


fa) =O (bd), 
Ig 


te} 4 s #0 (a). 


a. Proc. Paper 1706, July, 1958, by M.N. Fialkow. 
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(2) Along the simply supported ends y = + 1, 


(22) 
On page 1706-7, on the 4th and 6th lines, substitute ©” ‘in lieu of £ 


vq 
On page 1706-7, on the 24th line, in lieu of “Section IV”, substitute, “the 


section on Kinematically Admissible Velocity Fields”. 


On page 1706-12, delete the last three equations (63) and substitute, the 
following in lieu thereof: 


"Ga = G(o,0) = “By (0) -mo) -n,(o) < 1 
Ha = H(o) = -m(o) -n,(0) <1 (63) 


On page 1706-13, on the 2nd line, in lieu of “Section IV”, substitute, “the 
Section on Kinematically Admissible Velocity Fields”. 


On page 1706-15, on the 8th line, in lieu of “Section IV” substitute, “the 
section on Kinematically Admissible Velocity Fields”. 


On page 1706-15, on the 32nd line, substitute, 2™ “in lieu of cm ; 


On page 1706-18, on the 19th line, in lieu of “Section III”, substitute, “the 
section on Statically Admissible Stress Fields”. 


On page 1706-21, in the first of equations (58), the sign of the fourth term 
in the bracket in the numerator should be +. 


On page 1706-20, the 20th line should read “where A = 1 « 


On page 1706-22, on the 3rd line from the bottom, change the first word 
from “Considered” to “Certain”. 
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